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ABSTRACT 


The difference in magnitude between ¢ and « Persei has been determined throughout 
a spectral region of more than 5000 A with a photoelectric photometer used in con- 
junction with an objective grating. 

When these observations are combined with those by Kienle and by Mrs. Rudnick, 
the difference in magnitude between these two stars is, to a first approximation, pro- 
portional to the frequency. These data therefore do not definitely disprove the hypothe- 
sis that the peculiar color of ¢ is either partly or entirely due to a temperature differ- 
ence between these two stars. If this coloring is attributed to matter in interstellar 
space, the wave-length exponent in the attenuation term is near unity in the observed 
spectral region. A more literal interpretation of the data suggests that this exponent 
increases from about 1.2 in the visual to 1.6 in the infra-red. 

Infra-red colors and magnitudes of 19 B stars have been determined. A comparison 
of the color excesses computed from these data with those given by Stebbins and 
Huffer also indicates that this exponent is near unity. 


It has been shown in a previous paper’ that certain early-type 
stars measured in the infra-red region of the spectrum were much 
redder than the avertge star of the same spectral type. The same 
stars were known to be highly colored, from measures in the blue 
and yellow regions of the spectrum. 

In an effort to determine the attenuation coefficient as a function 
of the wave-length, the effective energy of ¢ Persei was compared 
with that of € Persei over a long spectral range, on five different 
nights at the Sproul Observatory. 

Since this program was started, Mrs. Rudnick, at the Yerkes Ob- 
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servatory, has published a paper? on this subject and has based her 
conclusions mainly on her spectrophotometric observations of these 
two stars. Kienle* has also made spectrophotometric comparisons 
of these stars. Three similar comparisons have been made in the blue 
and the red regions by the Greenwich observers Greaves, Davidson, 
and Martin.‘ The Sproul data were derived in an entirely different 
manner and cover a much longer spectral region than former 
comparisons. 

Spectra were produced by a fine-wire objective grating. A spec- 
trometer in the focal plane separated the two first-order spectra into 
definite energy units of known wave-length and reflected them to 
the photoelectric cell. This method has been described in a previous 
paper.® 

Since these two stars are only a little more than eight degrees apart 
in the sky, and since they were observed near the zenith, the correc- 
tions for differential extinction were negligibly small. 

Since I was interested in the difference in magnitude of these two 
stars as a function of the wave-length, it was important to know 
whether the scale of the photometer changed with the wave-length. 
From former laboratory tests I found that the observed response of 
the cell was linear if an incandescent lamp was used as a source. 
Since the two stars are about equally bright visually, it seemed more 
important to check the scale at 10,000 A, where ¢ Persei is about 0.6 
mag. brighter than ¢ and the observed electrometer deflections are 
small. The scale of the photometer was checked in this region in the 
following manner: the two first-order spectra of ¢ Persei were fo- 
cused on their respective slits, which were displaced from the central 
image at such a distance that their centers corresponded to 10,000 A. 
The slits were 2 mm or 980 A wide. The response was then measured 
with first one and then the other slit covered, and again with both 
slits open. This process was repeated with no star in the field, in 
order to determine the response due to background and other effects. 
After the stellar observations were corrected for background, it was 
found that the total response from the separate slits equaled the 
combined effective energy, within the error of observation. This 


2 [bid., 83, 394, 1936. 4 Tbid., p. 702. 
3 M.N., 88, 700, 1928. SAp. J., 84, 372, 1936. 
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process was completed four times, and the mean scale factor was 
1.00 + 0.02 (probable error). 

Figure 1 shows the photoelectric response as a function of the 
wave-length, according to the observations of December 7, 1936. 
Since the slit-widths were 2.00 mm wide, the effective wave-lengths 
were not exactly the same as the wave-length of the center of the 
slit, which is the abscissa of the curves. It is evident from the num- 
ber of observations that each spectral region was measured twice. 
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Fic. 1.—The photoelectric response in eyepiece divisions per second is plotted 
against the wave-length of the center of the slits. ¢ Persei is represented by dots, and 
e Persei by circles. If the positions of the slits were numbered from 1 to 13, beginning 
at the left, these data were observed in the order 1, 3, 5, 7, 9, II, 13, 12, 10, 8, 6, 4, 2,1. 
Readings with the standard lamp in position, and also readings due to background 
were made at each end of the stellar sequence. 


These observations are sufficiently numerous to enable one to deter- 
mine by a graphical process the wave-length which equally divides 
the effective energy for each setting of the slits. When this wave- 
length is known, a smoothed curve through the points may be used 
to determine the corrections to the observed energy in order to re- 
duce it to its value at the wave-length corresponding to the center 
of the slits. I have carried through this procedure and found that 
these corrections are 0.04 mag. or less, except for the last three points 
in the infra-red. Since we are primarily concerned with the difference 
in magnitude between these two stars as a function of the wave- 
length, the corrections to be applied to the measured magnitude dif- 
ferences are very small, except for observations beyond 10,000 A. 
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These differential corrections are given in the next to the last column 
of Table 1. The wave-length of the center of the slits, which in view 
of these corrections may be considered as the effective wave-length, 
is given in the first column of Table 1. 

The small grating mentioned in the footnote of Table 1 is the 
one mentioned in the former paper.s The large grating gives 


TABLE 1* 





MAGNITUDE DIFFERENCES: { PERSEI—e PERSEI 





WaAVE- 
LENGTH 
(um) 1936 1936 1936 1936 1936 |Weighted Weight Diff. Final 
Jan. 25 | Feb. 8 | Nov. 29 Dec. 7 | Dec. 22} Mean Corr. Mean 
EE ORR ec ae +-oM28 |.......|.... +oM28| 0.5 |...... 
©.477 +oMar |+0™16 + .31 |\+0™M2q4 +018 + .21 | 4.0 |—oMo2 |+o0™1g 
Ge Ree 6 ce es ae + .c6\— .o!I cre + .02 5.5 .C0 + .02 
O:573..0.+.)— 660 |— «Os IQ im .I2|— .O7 |— .10| 4.0 \- .OF |— .00 
et are .17 .14 eT ee ee | 2.5 3 -02 16 
0.669.... .24 .24 .28 me key 34 a5 | 4.0 00 .28 
oy Se Seed .25 <a5 Sd eae .28 25.5 .0O 28 
0.705.... -37 32 -37 -39 -42 37 | 4.0 oo 37 
ne ee Bese .48 Re a -43 2.5 .00 .43 
0.862 .50 .40 53 52 42 48 | 4:0|— .O% 49 
Se eee ONE ae ee-| i(.26) 1S eae .57| 2.0/— .O1 .58 
GOse: 25... 0.51 0.64 a2 .60 0.54 59 | 6.5'\r .O% .56 
O:002..... Se) Peres OSO0 Heist. BP 66 | 0.5 .0O 
I .007 mynhet ‘ae 52 ae 52 1.0 |\+ .o1? 0.62 
1.055. en 0.04 0.04 I.0 |—0.08 


Slit-widths. 2mm > 2mm/!1mm|{2mm_¥2mm 


Grating....| Small | Small | Large Large | Small 
Number of 
readings I 2 I 2 2 





* No observations were made on January 25 of the region beyond 0.8624. The value o.s51 at 0.9584 was 
derived from a recent comparison in this spectral region. The value in parentheses has been assigned zero 
weight in forming the mean magnitude difference. 

a similar dispersion at the focus but covers 0.9 instead of 0.4 of 
the objective. The hvdrogen emission lines Ha and H@8 in y Cas- 
siopeia have been seen and photographed with this grating in front 
of the 24-inch objective. The dispersion as derived from a single 
exposure of this star appears to be near 483 instead of 480 A per mm. 
The use of the former value, which was used for the smaller grating, 
does not noticeably change the results of this paper. Except for the 
observations of February 8, the focal length was changed in such a 
way that the slits were always in reasonably good focus for the par- 
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ticular wave-length bands passing through them. On February 8 the 
focus was fixed at a mean value for all readings. The two nights on 
which only one reading was made for each position of the slits were 
given half-weight in forming the mean magnitude differences. In 
general, since these single readings average longer than 30 seconds 
in duration, they are equivalent to at least two short readings. Ex- 
cept for these two nights, when the telescope was pointed at each 
star but once, the first star was always reobserved. The transparency 
was excellent on all five nights. 

Both ¢ and e Persei have companions which share their proper 
motions of 13” and 9” separation, respectively. Their respective 
visual magnitudes are 9.3 and 8.3. Both are on the edge of the field 
striking the cathode of the cell. Since, according to descriptions in 
Burnham’s catalogue, neither companion is reddish in color, it is 
highly improbable that they influence the results by an appreciable 
amount. 

The final mean magnitude differences of Table 1 are plotted 
against the frequency in Figure 2. If ¢ and ¢ radiate as black bodies 
with temperatures of 22,500° and 10,000° K, respectively, and if the 
theoretical curve were fitted to the straight line at its mid-point, the 
extreme values of this curve would fall 0.04 mag. below the straight 
line at each end. These data do not definitely refute the hypothesis 
that the peculiar color of ¢ is due to an unusually low temperature 
of this Br star. 

The color of € Persei has been found to be normal for its spectral 
class in both the violet and the infra-red regions of the spectrum. 
The spectral type and absolute magnitude found by Edwards? for 
¢ Persei are Biss and —3.2. The corresponding data for ¢ Persei are 
Bis and — 2.9. Rimmer’ has classified both ¢ and ¢ Persei as Bis and 
has assigned an absolute magnitude of — 2.5 to both stars. However, 
Struve and others have shown that ¢ Persei has strong interstellar 
lines, and is therefore very luminous. Plate III shows that the spec- 
tra of these stars are very similar. It is therefore reasonable to 
suppose that the distribution of energy is the same function of the 
wave-length for each of them, and that the observed reddening of 
¢ is due to selective absorption or scattering of its initial energy by 


6 M.N., 87, 377, 1927. 7 Commonwealth Solar Obs. Mem. 2, 1930. 
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particles in interstellar space. If we let e~** “t¢ be the transmission 
coefficient of intervening particles, then the observed difference in 
magnitude between these two stars is 


Am, = — (2.5 log e)kRN*7 + D. (1) 


The coefficient & is a function of the properties of the intervening 
medium and is assumed to be independent of \. According to 
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Fic. 2.—The observed magnitude differences are plotted against the frequency. 
Observations by Kienle are represented by small circles, and those by Mrs. Rudnick 
are denoted by dots. The large circles are the final mean differences given in Table 1. 
No systematic adjustment has been made between the three sets of data. Kienle has 
attributed the two discordant observations at the extreme left of the diagram to un- 
usually heavy absorption near the head of the Balmer series in e Persei. The data of 
the Greenwich observers are not shown. Their mean gradient of 0.62 agrees very closely 
with the value 0.61 indicated by the straight line in this figure. 


Schoenberg and Jung,* this assumption is valid for gases and non- 
absorbing spheres and is a much better approximation in the red 
than in the violet for small metallic particles. The constant D is 
independent of A, and its value depends on the luminosities and dis- 
tances of the two stars, as well as on the general absorption. 


8 A.N., 253, 266, 1934. 
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The data are well represented if a is set equal to unity. The value 
of kis then —o.61, and D equals —1.25 mag. If the slight upward 
curvature indicated in Figure 2 is real, a increases with \ and its 
rate of increase depends on the value assumed for D, as Mrs. Rud- 
nick has shown. I have recently discussed these data from this point 
of view in a paper read before the American Astronomical Society. 
One would not expect D to exceed the lowest observed Am, namely, 
—o.6 mag. In view of spectroscopic and other data, we may reason- 
ably assume that D is —1.0 mag. In this case a increases from 1.2 in 
the visual region to 1.6 in the infra-red. Schoenberg and Jung, as 
well as Schalén,? have shown by theoretical reasoning that a may 
have a value near unity. Gerasimovic’® was apparently the first to 
suspect this to be the case. Struve, Keenan, and Hynek"™ have found 
evidence that a is near 2. Schalén? finds that the \ law fits the ob- 
servational data in the photographic and visual regions better than 
Rayleigh scattering or the \~4 law. Since Am is near —o.6 mag. at 
10,000 A and +0.3 mag. at 4250 A, the total photographic absorp- 
tion is at least 0.9 mag. 

The unusual color of ¢ Persei has been discussed from two points 
of view: as a low-temperature phenomenon and as the result of se- 
lective absorption or scattering in interstellar space. Kreiken’s” re- 
cent suggestion that a considerable portion of the reddening in B 
stars may be attributed to stellar rotation does not seem to be valid 
in the case of this star. Miss Westgate" finds that the rotational 
velocity of ¢ is slightly greater than that of ¢. 

It seems important to separate stellar from interstellar reddening. 
It may be possible to do this by carefully studying the colors and 
spectra of the companions of B stars which share the proper motions 
of their primaries. In the case of ¢ Persei, for example, since the 
early B stars are among the bluest stars known, we would not expect 
to find its companion to be bluer than ¢, which appears yellow. If 
this were the case, it would imply that some of the reddening is of 
stellar origin. According to Burnham’s catalogue, this companion 
is “ashen.” It is very difficult, with my present photometric equip- 

9 Medd. Upsala Obs., 64, 1936. 

10 Zs. f. Ap., 4, 265, 1932. 2 A.N., 261, 202, 1936. 

1 Ap. J., 79, 1, 1934. 13 Ap. J., 77, 142, 1933. 
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ment, to measure the color of this companion. HD 170740 is another 
reddened B star with a ninth-magnitude companion. 

In 1934 Stebbins and Huffer’ published the results of an extensive 
color survey of early B-type stars. The effective wave-lengths of 
their colors are in the blue and the violet regions of the spectrum. In 
October, 1935, I observed the colors and magnitudes of 19 B-type 
stars which had been selected on the basis of their results. Eight 
comparison stars of later spectral types were also observed. The ef- 
fective wave-lengths of these colors are 6110 A and 8000 A, and of 
the magnitudes, 8000 A. These wave-lengths were derived from con- 
temporaneous observations of Vega with the objective grating. 
Since these B stars were all observed at small zenith distances, their 
colors and magnitudes were reduced to the zenith by means of the 
extinction coefficients formerly found for Ko stars.5 

The results of these observations are shown in Table 2. The vari- 
ous columns are, respectively: the HD number of the star, the mean 
color, the number of nights on which the star was observed, the 
visual magnitude, the infra-red magnitude, the spectral type, Steb- 
bins and Huffer’s color excess, and a color excess £, derived from the 
filter measures at effective wave-lengths 6110 A and 8000 A. A third 
color excess, E,, given in the final column is derived from colors 
formed by taking the difference between the visual and infra-red 
magnitudes. For the stars observed more than once the average de- 
viation from the mean for the observed colors is 0.02 mag., and for 
the magnitudes, 0.03 mag. 

Since Ko stars whose infra-red magnitudes have already been 
published were observed on the same nights as were the B-type stars, 
the zero point and magnitude scale of these stars are the same as for 
the Ko stars. In forming the color excesses, the mean colors of each 
spectral type were derived from a curve involving the comparison 
stars and B stars whose color excesses, as given by Stebbins and 
Huffer, are less than +0.04 mag. The mean differences in magnitude 
were derived in the same way with the same stars in the deriva- 
tion of Ep. 

The two observed magnitudes of HD 223960 are very discordant. 
Evidently on one of these nights a near-by K5 star (HD 223866) of 
similar visual magnitude was mistaken for the Bo star. As I had 


14 Pub. Washburn Obs., 15, Part 5, 1934. 






















































B-STARS IN THE RED AND INFRA-RED 153 


relied entirely on the circles of the telescope in locating these stars, 
I have recently looked up the visual magnitudes and spectra of 
the stars in the vicinity of those observed B stars which are fainter 


TABLE 2 















































MAGNITUDE Cotor Excess 
STAR , ee er a ey ee 
- q es Infra- ” 
Vis. Red Ry E. Em 
3940... .....| +018 I 7.40 | 6.04 Bs | +o™go | +0™63 | +1™q43 
PO ion. 55d tte ee I 6.93 | 6.62 Bo 14 20 | +0.59 
7 eer 38 I 6.08 | 6.07 B7 fore) .04 | —0.01 
14818. . i I 6.2 5.67 Bi 27 36 | +0.81 
206165... ay 25 2 4.87 | 4.2 B2 .22 25 0.82 
206672. 48 2 4.78 | 4.68 B3 03 oo 0.25 
QO7 EOF. i. 48 I 6.16 | 6.15 B3 03 .0O 0.16 
207708......: 18 I 6.56 | 5.90 B2 . 26 32 0.85 
210830. . 25 I 5-19 | 4.59 06 19 | + .31 | +1.04 
212222.. 49 I 6.27 | 6.34 B7 eco | — .07 | —0.09 
214993. 50 3 g-¥8 | 44 Br oo | + .o1 | +0.10 
216014. i 3 6.83 | 6.34 Bo 27 35 0.77 
215200.. 34 2 5.84 | 5.56 B3 12} + .16 0.43 
216010 :.......... 49 I 5. s4. [5:61 B3 o2 | — .OI 0.08 
2IG349.... : .06 I 7.40 | 6.75 B2 30 | + .44 | +0.90 
73) ee 49 I ey Oy ee B3 OI | — .OI | —0.09 
36. Se .45 I 4.93 | 4.81 Br 03 | + .06 | +0.36 
SERA IO si 6-06 e165 — 66 2 5.84 | 6.06 B3 03 | — .02 | —0.07 
QAZ0000 . is sc cel oe Ed I 6.98 | 5.69 Bo | +0.43 | +(0.63)| +(1.57) 
223866? .......-.. +o.18 I 6.92 | 5.30 eee bas 5 fcis Bee pe) fe sh ee 
Comparison Stars 
| 
872: —0O.15 I 5-32 | 4-75 | Fr | | 
8763... | + .14 I 5.63 | 4.25 Ki Ae Soavacy eae 
11443... ..| — .18 I 3.58 | 2.78 | F2 | 
11973.. . 26 I 4.83 | 4.33 | AS | AE Cay 
12869... a 28 I 5.08 | 4.95 | Ao ER Tee 
0? ee 38 I 4.85 | 4.78 | Az }-- cperarmegee| Se Fe? 
215182.. eee .00 I 3-10 | 2.25.) Ge: f.. hive ss 0dabedaeuaen 
2900RF 2s 66 ..| —0.09 I 4.57 | 3 76 | F6 | oe 











than 6.00 mag. It appears highly improbable that any other B-star 
was misidentified. If we let C be the observed color and M, the ob- 
served infra-red magnitude, then the visual magnitude may be com- 
puted by the empirical formula 


M, = M, + 2.30C + 1™o2. 
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The average difference between the computed and the observed 
visual magnitudes for the B stars is 0.09, and the largest difference 
is 0.18 mag. This procedure affords a partial check of the identity of 
the stars, since the visual magnitudes were taken from the HD cata- 
logue. 

In the notation of equation (1) the color excess 


E = Am, —Am, = — (2.5 log e)k(Ar* — AZ), 


where A, >A,. If we let E, be the color excess of a star as observed 
by Stebbins and Huffer, the ratios 
: taeda Ee A 


= <a nes 
i a a En A7~* — ATe’ 
d 5 3 ] 


where A, = 0.426 w, A. = 0.477 M, A; = 0.550 mw, A, = 0.611 w, and 
A; = 0.800 uw. Since the effective wave-lengths are known to a first 
approximation, we can compute these ratios and compare them with 
the observed values for different values of a. The mean values of E,, 
E.,and E,, have been used in the computation of the observed ratios. 


The results are given in Table 3. 





TABLE 3 
CompuTeD RATIO 
ae eee 
| E;/Ec ee ae 
4.0 2.34 I.30 
2.0 I.00 0.64 
1.0 0.05 0.44 
0.6. 0.55 0.38 
Observed 

ratio 0.71 0.28 





In making this comparison it is well to remember that we have 
again made the assumption that stars of the same spectral type have 
the same energy distribution and that & is not a function of A. We 
make an additional assumption that the effective wave-lengths are 
independent of the color. This can only be regarded as a first ap- 
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proximation. The effective wave-lengths used for the color excesses 
found by Stebbins and Huffer" have not been corrected for two com- 
pensating factors: the energy distribution of early-type stars and 
the atmospheric extinction. 

The data in Table 3 indicate that the computed ratios of color 
excess agree with the observed ratios when a is near unity. £, is 
probably more reliable than £,,. Within the errors of measurement, 
E,/E. is constant for all reddened B stars in Table 2. Although we 
may not know the true cause of reddening, the fact that the data in 
Table 2 are consistent within themselves and agree well with those 
found for ¢ Persei suggests that similar physical processes are in- 
volved for a number of stars. 

SPROUL OBSERVATORY 
SWARTHMORE COLLEGE 
SWARTHMORE, Pa. 
February 1937 








LINE CONTOURS OF THE ATMOSPHERIC 
OXYGEN BANDS* 


Cc. W. ALLEN’ 


ABSTRACT 


The breadth and structure of the atmospheric oxygen lines have been determined 
by studying the curves of growth of the three bands A, B, and a in the solar spectrum. 
True central intensities are thus computed as a function of the equivalent breadth. The 
absolute total absorption coefficients for the three bands have also been found. 


INTRODUCTION 

A considerable amount of experimental photometry of the atmos- 
pheric oxygen lines has led eventually to the conclusion that they 
are broadened essentially by molecular collisions and that the ex- 
ponential law of absorption in the atmosphere is obeyed.” * #5 
Hence the complete contour of a line of a given intensity may be 
computed if the main constant of collision broadening (e.g., the 
effective collision diameter of the molecule) is known. Since such 
computed contours are narrower than solar lines but similar in other 
respects, they are ideal for determining the performance of a spectro- 
graph when the contours of Fraunhofer lines are to be measured. 
The main broadening constant, however, has been determined only 
indirectly,’ and hence a redetermination is necessary before the 
lines can be used for photometric control. The method used in the 
present paper is that of studying the curve of growth, i.e., the curve 
relating the logarithm of the equivalent breadth to the logarithm 
of the number of absorbing atoms. 

THEORY 

Theoretical curves of growth have been computed by using fixed 
collision broadening factors. The number of collisions per second, and 

* Contributions from the Mount Wilson Observatory, Carnegie Institution of Wash 
ington, No. 566. 

* Hackett Research Student. 

2?R. v.d. R. Woolley, Mt. Wilson Contr., No. 420; Ap. J., 73, 185, 1931. 

3W. H. J. Childs, Ap. J., 77, 212, 1933. 

4 Ibid., Phil. Mag., 14, 1049, 1932. 

5 H. von Kliiber, Zs. f. Ap., 6, 161, 1933. 
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hence the collision broadening, is proportional to the atmospheric 
pressure (the effect of variation in temperature being comparatively 
small). Suppose the collision breadth of the absorption coefficient 
to be 2b(p/p.), where 2) is the breadth at normal pressure, and p 
and p, are the pressures at height / and at ground level. The ab- 
sorption per centimeter for a line is 


a b (p/ po) _?p 
x (v — vo)? + 8 p?/ps Po 


ay = ’ 
where a is the total absorption for the line at normal pressure. Put 
pb = pe”, where H is the height of the homogeneous atmosphere; 
then 

a be—2/H# 


- (v am Yo)? + be 2h/H ? 


and the total absorption through the whole atmosphere is 


Ae 


ai 
ab e2"/H sec z dh 
T (v — )? + be 2h/H ” 
(1) 
b? 


= aH sec z : 
7 (? — Wa)" | 


ae — log. (1 + 
where. z is the sun’s zenith distance. 

This contour has a peak at vy — v, = 0, but it is smoothed out in 
practice by Doppler broadening. The contour (1) was combined 
with Doppler broadening .(assumed temperature 270° K) by nu- 
merical integration. With the absorption coefficients from the con- 
tour so obtained, absorption lines were plotted and the areas meas- 
ured to derive equivalent breadths. The process was repeated for 
different values of a to obtain the curve of growth. 

In order to derive the experimental curve of growth, it is necessary 
to determine the relative intensity of various lines in the band, i.e., 
the relative values of a. This is given by 


a = CNie—E/t? , (2) 


where C is a constant, NV the number of molecules capable of ab- 
sorbing the band per cubic centimeter, e~“*” the Boltzmann factor, 
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and 7 an intensity factor which is known for the oxygen bands.° 
The temperature of the atmosphere is irregular, but it is sufficient 
for our purpose to adopt the formula 

5-3 SC 10° 


TT. 
h+ 5.3 X 10° 





where / is the height in centimeters and 7, the ground temperature. 
The absorption for a line at any height is then given by 
eer 
a, = CN,ie~E/*T — , 
¢ 0 
where J, is the value of V at normal temperature and pressure, and 
d and d,, the densities of the atmosphere, are related by 


d 
dy 


=_ enh QX 105 


The mean value of a to be used with a homogeneous atmosphere is 
then 


5 = CNet { Re de mee _ a ee 
~ 9 X 105 J, sil kT, s.3 * 10° 9 X 105} 





which, except for the denominator, has the same form as (2). A 
somewhat similar expression was used by Mecke and Childs.’ 
OBSERVATIONS 
A number of plates were secured with the tower telescope and 
75-foot grating spectrograph of the Solar Laboratory of the Mount 
Wilson Observatory, situated in Pasadena at an altitude of about 
700 feet. Photometric standardization was attained by use of step 
weakeners which had been calibrated at the Physical Laboratory 
at Utrecht. The weakeners were kindly lent by Dr. Mulders. 
6W.H. J. Childs and R. Mecke, Zs. f. Phys., 68, 344, 1931. 
7 Ibid., p. 362. 
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TABLE 1 


INTENSITY DATA FOR THE OXYGEN BANDS A, B, AND a 
































EQUIVALENT BREADTH IN cm * 
, p 
LINE eae a eek ei “ Re log a 
cm 
A B a 
P; (0). stachinge I.10 0.50 ©.100 3 0.00 
P;.¢ 2). ; 1.24 56 | 126 15 0.06 
P; ( 2) we valte ede 66 | 150 17 0. 26 
P;( 4) 1.50 | 62 | 156 41 0.25 
P, ( 4) | 1.65 | 77 178 43 0.37 
P; ( 6) 1.59 67 166 79 0.33 
P; ( 6) roa 70 188 81 0.41 
P; ( 8) : 68 168 127 0.32 
P; ( 8) ; 72 | IQ! 129 0.39 
P; (10) ; ‘ 63 } 156 188 0.27 
P, (10) : | .67 .164 190 0.33 
P; (12) | 55 133 260 0.18 
P, (12) I. 39 58 | 140 262 ©.22 
P; (14) SS mae 49 | 121 343 0.05 
Pa trays. See pase ee 52 | 121 345 0.09 
P (16) “k .F:66 4o. | 085 439 —o.10 
P (18) 33. | ~~ ~-060 545 —0.30 
P (20) 0.59 ate 040 662 —0.52 
P (22) fn 0.46 .170 c20 791 —o.78 
P (24). Be ee 128 0126 931 —1.06 
P (26) : 0.21 061 | 0.0075 1083 —1.37 
P (28). acel Oka 038 | er 1246 —1.71 
P (30) eer 0.104 0.023 | ; 1420 —2.05 
P (32) os ed 0.042 1606 —2.46 
| 
P (34) 3 zt 0.0456 1803 — 2.86 
POGGh oes atu eee 0.0064 | 2011 — 3.30 
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Table 1 gives the observed values of equivalent breadth in wave- 
number units for the P, and P; branches of the bands A, B, and a 
situated at 7600, 6900, and 6300 A. The nomenclature is that used 
by Badger and Mecke.* Since up to J = 14 there is a measurable 
difference between the intensities of the two branches, both values 
are listed. The table also gives the excitation potentials in wave- 





log equivalent breadth (cm-) 

















log a 


Fic. 1.—Curves of growth for atmospheric oxygen lines 


number units,’ and relative values of @ computed from (3) and based 
on a ground temperature of 303° K. The actual mean temperature 
during the period of observation was 301° K, but this was increased 
slightly to allow for the peculiarities of the local temperature distri- 
bution. The air masses in terms of a homogeneous atmosphere at 
sea level were, respectively, 1.2, 2.9, and 5.0 for the bands A, B, and 
a. The observed curves of growth are plotted as points in Figure 1; 
8 Zs. f. Phys., 60, 59, 1930. 


9 Obtained from the measurements of G. H. Dieke and H. D. Babcock, Proc. Nat. 
Acad., 13, 670, 1927. 
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and, when fitted to a set of computed curves, the following values of 
b were obtained: 


Value of 6 (cm~') 0.049 | ©.050 








The mean value of 0.049 cm™ has been used for computing all the 
curves shown in Figure r. 

Two further checks have been made on this value of 6. The a 
band has been measured in a solar spectrum taken on Mount Wilson 
(altitude 5715 feet); and, after correction for the decreased pressure, 
it gave a value of 0.051 cm~'. The A band has also been measured 
in a 475-foot air-path spectrum taken on Mount Wilson with a 
tungsten strip lamp as source of light. In this case the pressure and 
temperature were practically constant over the length of the air 
path, so that the curve of growth is simpler to compute. It is very 
satisfactory to find a value of 0.048 cm~ for 6. The agreement with 
the values found for the whole atmosphere gives strong support to 
the assumption made above that the broadening increases with 
pressure. 

CONTOURS OF ATMOSPHERIC LINES 

When the value 6 = 0.049 cm is introduced into (1) and a 
Doppler broadening equivalent to T = 270° K and molecular 
weight 32 is superimposed upon this contour by numerical integra- 
tion, the resulting frequency distribution of A, is that given in 
Table 2. Actual contours of exponential absorption lines with this 
value of A, have been plotted, and the relation between equivalent 
breadth and central intensity so obtained is shown in Table 3. This 
table should be useful in determining the performance of a high- 
dispersion spectrograph when used with absorption lines. If the alti- 
tude of the observing station is considerable, the lines will be some- 
what narrower and the central intensities less. To obtain the value 
of the central intensity in this case from Table 3 the measured value 
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of the equivalent breadth should be divided by the reduction in 
pressure. 
PHOTOMETRIC DATA FOR OXYGEN BANDS 
The observations enable one to determine the following photo- 


metric constants of the oxygen bands: (a) the effective collision 
diameter of O, in air; (b) the relative absorption coefficient of the 


TABLE 2 TABLE 3 
CONTOUR OF THE ABSORPTION Co- RELATION BETWEEN CENTRAL INTEN- 
EFFICIENT Av OF OXYGEN SITY r, AND EQUIVALENT BREADTH 
ATMOSPHERIC BANDS W FOR OXYGEN BANDS 
v—ve (Cm—") Ay W(cm-") fe 

2 re AE SO A dd ciate  1OF4 
0.01 rr25 INE cr Av ae shesk wes eth 0.53 
6:02: EY j OrO0% 25 0.36 
0.03 Sad CIOS a6 os 0.22 
0.04 3.6 0.10 0.13 
0.05.. 2.4 0.12 0.07 
0.06. 1.8 0.14 0.04 
0.08 rr 0.16. 0.02 
0.10 0.70 o.18 0.01 
See ue ew 0.34 OO Owe ok eos 0.00 
0.20 0.188 

0.30 0.088 

0.40 0.049 

0.50 0.032 

I.00 0.008 

2.00 0.002 


three bands A, B, and a; and (c) the absolute value of the absorption 
coefficient. Most of these constants have been studied already by 
Childs* 4 and hence will receive only brief mention. 

a) From Dennison’s formula for broadening due to collisions 


10 


we 
may write 


b(sec-') = 1.10N0? Vaz : 


where 7 is the number of molecules per cubic centimeter, and o the 
molecular diameter, and where k, 7, and M have their usual mean- 


10 Phys. Rev., 31, 503, 1928, formula (1”). 
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ings. Introducing the numerical values, m = 2.7 X 107, k = 8.32 
X 107, T = 290 K, M = 30, b = 1.46 X 10° sec, we find 


o = 4.2 X10°cm. 
This result may be compared with Childs’s value obtained indirectly: 
o=2.5X 10 °cm. 


b) The relative absorption coefficient in the three bands, when 
corrected for the difference in air mass and expressed in wave- 
number units, is 


This result is in agreement with the work of Childs, whose corre- 
sponding factor for the ratio A:B is 11.2. 

c) Figure 1 includes a curve relating the equivalent breadth to 
the absolute absorption coefficient integrated over the breadth of 


TABLE 4 








Absorption 


Coefficient Molecular 
per Centimeter Absorption 
of Normal Coefficient 


Atmosphere 


cm—'sec™? cm?sec! 


7620.987 A... 4.2X105 | 7.3Xxk10°% 

Total A band... 2.7X107 4.7K10> 

Total B band... 2.0X10° | 3.5X107% 

Total a band. | 7.8X104 1.4X1074 
| 





a line in wave-number units. Thus, for any measured line the ab- 
sorption coefficient for the whole atmosphere may be read from this 
curve. Dividing this coefficient by Hsec z and multiplying by c, 
the velocity of light, we find the absolute absorption coefficient in 
frequency units for an atmosphere of normal temperature and pres- 
sure. The molecular absorption coefficient is obtained by dividing 
again by the number of O, molecules per cubic centimeter at normal 
temperature and pressure. The results are as shown in Table 4. 
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The total absorption per centimeter at normal temperature and 
pressure for the A band has also been found from the 475-foot air- 
path plate. The result is 2.6 X 107 cm™ sec™, in excellent agree- 
ment with the value obtained from the whole atmosphere. Childs 
and Mecke’® obtained the value 4.8 X 107 cm~ sec~ for this con- 


stant. 


In a forthcoming paper the oxygen-line contours will be used to 
determine the central intensities of Fraunhofer lines. 
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CENTRAL INTENSITIES OF FRAUNHOFER LINES* 
C. W. ALLEN’ 


ABSTRACT 

The central intensities of gt Fraunhofer lines in the red and infrared regions of the 
solar spectrum have been measured. The corrections for resolving power were obtained 
from a study of the atmospheric oxygen lines which appeared on the same plates as the 
solar lines. 

By plotting the central intensity against the equivalent breadth it is found that all 
observed lines are consistently wider than theory predicts. In the main, the discrepancy 
may be overcome by assuming that the sun’s reversing layer is subject to turbulent 
velocities of the order of 1.6 km/sec. The discussion deals with the nature of the turbu- 
lence, the change in the structure of Fraunhofer lines with wave-length, and the causes 
of abnormal sharpness and breadth in individual lines, and leads to the conclusion that 
the observations may be explained without departing from the simple method common- 
ly used in computing Fraunhofer contours. 

The intensity features of a Fraunhofer line may be described 
conveniently by the equivalent breadth and the central intensity. 
Numerous photographic observations of equivalent breadth have 
been published; but in order to observe central intensities by photo- 
graphic photometry, one must overcome the difficulty of obtaining 
a suitable correction for the distortion of the absorption line by the 
spectrograph. Carroll? and Plaskett’ have computed the corrections 
for ideal spectrographs giving theoretical resolving power, no ghosts, 
and no scattered light. Redman‘ and Thackeray’ have eliminated 
ghosts and scattered light by using a monochromator. Shane® has 
employed interferometry to determine the correction, while the 
present author7* has used fine emission lines. It has not been pos- 
sible, however, to test the corrections so found on lines whose con- 
tours are known, and the disagreement between former results 
shows that the errors are considerable. In the present paper the 
corrections to be applied to solar lines are derived entirely from a 

* Contributions from the Mount Wilson Observatory, Carnegie Institution of Wash- 


ington, No. 567. 


* Hackett Research Student. 4 Ibid., 95, 290, 1935. 
2M.N., 88, 154, 1927. 5 [bid., p. 293. 
3 [bid., 91, 870, 1931. 6 Lick Obs. Bull., 16 (No. 449), 76, 1032. 


7 Mem. Commonwealth Solar Obs., Canberra, No. 5, Part 1, p. 10, 1934. 


5 M.N., 96, 508, 1936. 
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study of atmospheric oxygen lines whose contours have recently 
been determined by the author;? and as a result, the effects of re- 
solving power, light-scatter, and even Eberhard effect are eliminated. 


OBSERVATIONS 

Spectra of the center of the sun’s disk have been obtained with 
the two 75-foot plane-grating spectrographs of the Mount Wilson 
Observatory. A Michelson 8-inch grating was used at the Solar 
Laboratory in Pasadena, and the Michelson 66 grating on Mount 
Wilson. Both spectrographs are fed by tower telescopes with equiva- 
lent focal lengths of 150 feet. The first, second, and third orders 
have been used in various regions, with dispersions ranging from 1.3 
to 5.3 mm per angstrom. The theoretical resolving power is about 
170,000 in the second order for both spectrographs. Eastman III F, 
144 S, IN, and 144 N backed plates were used and were brushed 
during development with metol-borax for 5 minutes at 18° C. The 
developer was diluted considerably in order to decrease contrast, 
chemical fog, and Eberhard effect. Step-wedges, kindly lent by Dr. 
Mulders, were placed on the slit to give the plate-calibration. These 
wedges had been calibrated at the Physical Laboratory at Utrecht, 
but were recalibrated on Mount Wilson by photographic comparison 
with a step-slit; the two calibrations were in good agreement. The 
regions studied were 6122-6336, 6828-6925, and 7657-7788 A, the 
choice being such that atmospheric oxygen lines, as well as solar 
lines, appeared on the plates. The spectra were measured photo- 
metrically in the usual way with the Mount Wilson microphotometer 
constructed by Dr. Dunham. For each line, both solar and atmos- 
pheric, values were obtained of 7,,, the uncorrected central intensity 
(with continuous background as unit), and of A,,, the measured 
breadth of the line at half its maximum depth. Using the results of 
an earlier paper,? we may obtain r,, the true central intensity, for an 
atmospheric line from the equivalent breadth W, which for prac- 


tical purposes is given by W = 1.3 An (1 —7n). The ratio of the 


measured to the true depth, i.e. (1 — r,)/(1 — r-), represents the 
instrumental correction for atmospheric lines. The value of this 


9 Mt. W. Contr., No. 566; Ap. J., 85, 156, 1937. 





































CENTRAL INTENSITIES OF FRAUNHOFER LINES 167 


ratio, plotted as a function of A,,, defines a curve which may be used 
to determine a first approximation for the correction to solar lines. 

The true central intensities of all atmospheric lines that are strong 
enough to be of importance are zero, and hence this method is not 
dependent on a very accurate determination of the structure found 
for the oxygen lines. There may be a small difference, however, 
between the correction-curves for atmospheric lines and for solar 
lines, owing to the differences in the line contours. In order to study 
this question, it has been necessary to obtain some idea of the in 
strumental contour, which we denote by K(/), where / represents dis- 
tance on the plate. The central part of K(/) may be obtained from 
the computed resolving power and slit width of the spectrograph, 
but an appreciable wing intensity must be added to this result. By 
combining the estimated value of A(/) with the known contour of 
atmospheric lines /(/) and numerically integrating the standard 
equation 


J(l) = {Ko —1)I(d)d, 


. 


we obtain J(/), the observed contour; and by repeating the inte- 
gration for lines of different intensity, we may compute a theoretical 
corrective-curve. Adjustments are then made in K(J) until the 
theoretical curve agrees approximately with the observed curve. 
The final K(/) may then be combined with lines of solar type to ob- 
tain another theoretical correction-curve; and the difference be- 
tween these two theoretical curves indicates the adjustment that 
must be made in the observed correction-curve when it is to be used 
with solar lines. The adjustment has been applied, but our imperfect 
a priori knowledge of Fraunhofer contours leads to an ambiguity 
of 0.01, or perhaps 0.02, in the final values of 7. 

An indication of the statistical errors may be obtained from the 
variation in the individual values of r. and the scatter in the points 
from which the correction-curves were obtained. The solar lines 
have been measured about ten times. In the 6300 A region the O, 
lines were not strong enough to provide values of A,, as large as 
those of the solar lines, and the correction-curve had to be extra- 
polated by comparison with the curve at 6900 A. In view of this 








*% 
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extrapolation and the ambiguity already mentioned, it would not 
be safe to compute the probable error of r, from statistical considera- 
tions. This probable error is not likely, however, to be greater than 
0.02, which is a considerable improvement on earlier observations 
in these regions. 

The results obtained by the present method lie between the two 
sets of values already published by the author.”* Evidently the 
values obtained from the prism spectrograph at Canberra’ were 
undercorrected, a result which is not surprising; but the overcorrec- 
tion of the grating observations made at Cambridge® requires some 
explanation. In the first place, the contour there assumed for the 
solar lines is not entirely suitable, and the blunter cores that Fraun- 
hofer lines probably possess would necessitate a smaller correction. 
Secondly, the neon lines used to obtain the instrumental contour 
may not have been perfectly sharp. Neon spectra taken recently 
at Mount Wilson exhibit broader lines than can be attributed to 
instrumental effects. Whether this is due entirely to some feature 
of the particular discharge tube used at Mount Wilson is not cer- 
tain, but at least it throws some suspicion on the use of the neon 
spectrum for the determination of the instrumental contour. 


RESULTS 

The results are set out in Tables 1, 2, and 3, in which the first 
three columns are taken from the Revised Rowland. The following 
columns give values of r. obtained with different spectrographs: 
(a) first order and (b) second order, Michelson 8-inch grating, Solar 
Laboratory; (c) second order and (d) third order, Michelson 66 
grating, Mount Wilson; (¢) mean of earlier observations made at 
Canberra (prisms) and Cambridge (grating). Column (e) shows 
fair agreement with the other columns but has not been included in 
the final mean of r.. 

The values of the equivalent breadths W (unit, 0.001 A) have 
been averaged from all available material, including numerous re- 
measurements made in connection with the present research. For 
this purpose the formula W = pA,,(1 — rm) has been found very 
useful, where p is a factor which depends to a small extent on the 
spectrograph and the nature of the line. For most cases p is between 
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TABLE 1 


CENTRAL INTENSITY, EQUIVALENT BREADTH, AND SHARPNESS, 6122-6336 A 












































re 
ry Ele E.P a W 10° rn d 
ment | | mA 
(b) (c) | (d) e) | Mean 
6122.231...|Ca@ 1.875] 0.10 | 0.106 | | O.12 10.163 215 | 35.0 
6126.230...|Ti 1.062 78 77 | .98| .777 20 Ce re 
6127.918...|Fe ee 51 52°) . SEA 47 7.7 1 1.00 
6128.990...|Ni 1.669] .73 72 | .76 | 720 23 3.8 
6130.147...|Ni 4.248] .78 78 | 90 776 46) a Fr 
6136.631...|Fe 2.443 21 20 12 200 142 | 23.1 1.02 
6137.709...|Fe 2.577 20 22 14 205 138 | 22.6 | 02 
6141.733...|BatFe| 0.701 17 18 24 175 22 | 19.9 | I.10 
6142.499...|St 5.595 74 | | .736 35 | 5-7 | 0.72 
6149.255...|Fet 3.873 66 64 | .66 | .642 38 6.2 | 0.89 
6151.630...|Fe 2. 167 47 17 | 50 467 47 +61 2.86 
6154.235...|Na 2.003 70 | 701 35 5.7 | o.86 
6155.148...|.Sz aC | ae 48 | 483 84 | 13.6 | 0.76 
6157.739...|Fe Se Fa eee 41 | 41 410 60 | 9.8] 1.04 
6159.387...]|.. 87 872 12 | ee 
6160.759...|Na 2.095 55 | | 54 552 54 8.8 | 0.85 
6162.185...|Ca 1.891 | 13 | it 132 296. | 43: Othe oni: 
6165.369...|Fe Se | 54 532 45 7.3 | ¥.04 
6166.446...|Ca 2.510 .40 34 308 71 | 11.6 | 0.98 
6169.050...|Ca et 32 32 317 90 | 14.6 | 0.98 
| 
6169.570...|Ca 2.515 26 28 20 278 111 | 18.0 | 0.98 
6173.348...|Fe 2.213 35 36.10.38 | .35 1 .36 67 | 10.8 | 1.06 
6176.822...|Ni 4.071 AS 45 44 45 449 63 | 10.2 | 0.94 
6180.216...|Fe 2.716 45 45 46 47 454 53 8.6 | 1.04 
6186.725...|Nz 4.088 72 ee ee wee 27 ee ere 
6188.002...|Fe 3-926] .54 55 52 54. | -542 45 | 7-3 | 1.00 
6191.192...|Ni 1.669 35 36 35 34 355 7E | 11.5 | 3.00 
6191.577..|Fe 2.4321 .22 21 22 16 216 537 {22.2 | ¥2Ga 
6200.327...|Fe 2.597 36 <a0 | 34 48.) 464 68 | 11.0 | 1.02 
6213.443...|Fe 2.233) .31 29 .32 34 307 81 | 13.0] 1.04 
6216.366...|V 0.272 7 68 .70 72 692 1 5.5 86 
6219.294...|Fe 2.188} .28 27 28 30 | .277 88 | 14.2 |] 1.04 
6223.996...|Ni 4.088] .73 72 ay 72 722 26 ree 1 eee 
6226.745...|Fe 3.0071 .91 7° 73 724 940 26 ae ae 
6229.240...|Fe 2.833] .60 60 63 64 | .607 34 5.5 | 1.06 
6230.111...|Nz 4.088} .81 81 83 81 817 17 ak, it Ree 
6230.742...JFe V | 2.548) .23 aI 22 .18 | .218 166 | 26.6 | 0.098 
6232.655...|Fe 3.638) 35] -35| -37 34 | -354 83 | 13.4 | 0.97 
6339 3266 be ues a 59 59 62 4) . 807 61 9.8 | 0.73 
6240.659...|Fe 2.213] 0.50 | 0.49 | 0.49 | 0.48 |0.493 46 7.4 | 1.08 
| 
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TABLE 1—Continued 
Ve 
d wae ee W loon 
ment mA 
(b) (c) (d) (e) Mean 

6246.333...|Fe 3.587) 0.28 | 0.28 | 0.26 | 0.24 [0.273 124 | 19.8 
6247.569...|Fet SONS! 252 51 .50 .48 | .506 53 8.5 
6251.845...|V 0. 285 88 87 86 | .874 13 2.1 
6252.572...|Fe 2 ABA «<5 a0 2 23 .18 235 $23 1 1007 
B2EE O06. ILE fess or 78 80 .79 787 18 2.9 
6259.598...|Nz Pe | ee 84 84 .84 | .843 15 2.4 
6261.112...|7?# I.424].. 50 54 Jef | 0S¥8 48 7B 
6265.148...|Fe ee | | 30 .29 .30' | ..206 S21 23.3 
6270.237...|Fe 2.846 .44 .46 45 | .449 53 8.5 
6271.289...|Fe ec) ae 76 76 | 0.78 757 21 2:2 
6297.808...|Fe 2.213 30 34 35 331 orl £2.32 
6301.517...|Fe 3.638 20 28 ae .282 122 | 19.4 
6302.508...|Fe 2 G71)  <37 35 a .3§2 81 | 12.8 
6311.511...|Fe 2.819 70 70 73 713 2 3.8 
6314.675...|Nz 1.927 39 40 40 396 aL Vede 2 
6318.036...|Fe 2.443 30 29 28 2g1 104 | 16.4 
6322.701...|Fe 2.597 37 36 35 361 76 | 12.0 
6327.613...|Ni 1.669 60 62 ere 619 LG £.5 
6330.103...|Cr 0.937 68 70 7h a See 699 27 ae 
6320.859...)f@ }...... .66 66 5 al Ae 662 33 Le 
6335-345...|Fe 2.188 27 27 PBs san aus 372 100 | 15.8 
6336.837...|Fe 3.671] 0.30 | 0.29 | 0.31 0.298 108 | 17.0 
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TABLE 2 


CENTRAL INTENSITY, EQUIVALENT BREADTH, AND SHARPNESS, 





6828-6925 A 











Te 

(b) (c) (d) Mean 
6.50 | 0.53 | 0.52 {O.517 
70 13 I61 .7tO 
-45 .49 49 | .480 
.63 .66 .66 651 
7 49 .49 | .486 

72 -73 -73 725 
87 84 85 856 
.42 2 2 22 
0.52 52 -53 | -§22 
saree 0.69 | 0.70 |0.694 
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1.2 and 1.3; but for lines with very extensive wings, such as the 
solar hydrogen lines, it may be greater than 1.6. 

Tables 1, 2, and 3 also give the quantity W/X, which for some 
theoretical studies is more convenient than W itself, and, in the 
last column, d, which represents the difference in sharpness of 
Fraunhofer lines and might be called the “relative sharpness.” It is 


TABLE 3 


CENTRAL INTENSITY, EQUIVALENT BREADTH, AND SHARPNESS, 7657-7788 A 

















‘A 
N vend E.P. ‘s 10°W /d d 
ment mA 
(a) (c) Mean 

ORT OOS 6s 6-0:0 5 Mg 5.086 Ce a 0.522 100 13.1 0.85 
9680.207. <5 Mn 5.407 57 0.54 565 86 11.2 0.84 
SOGE. COP. 6 ssc 6: Mg 5.728 58 56 572 110 14.3 0.72 
7698.977... K 0.000 18 18 180 150 19.5 ¥. 22 
77%0.367...... Fe 4.202 55 54 551 65 8.4] 1.03 
ri 8 fy BY; Fet 3.887 67 65 663 47 6.1 0.95 
7714.309 Ni I.927 2 43 426 103 53.3 I.02 
7715-577 Ni 3.683 64 64 640 48 6.2 I.00 
7723.21 Fe 2.269 66 66 657 38 4.9 1.13 
7727.616 Ni 3.663 48 46 473 go 11.6 1.01 
7742.722.. Fe 4.970 56 57 563 103 13.3 0.77 
7748.284.. Fe 2.936 42 41 418 102 53.4 I.04 
7748.894. . Ni 3.690 47 47 469 87 41.2 1.03 
7751.117.. Fe 4.970 68 67 679 44 5-7 | 0.95 
7771.054..- O 9.106 .66 67 660 79 10.2 0.68 
7774.177 O 9.106 70 7° 696 71 g.1 0.66 
7775-394 O 9.106 75 76 75° 58 7-5 | 0.62 
7780.507.. Fe 4.454 42 2 423 124 10.0 0.93 
7788.933.. Vi 1.942 0.50 ©.45 | 0.490 go 11.6 0.98 
































defined as the ratio of the line depth 1 — 7, to the mean depth of other 
Fraunhofer lines of the same equivalent breadth and wave-length. 
These values of d are more accurate than those found hitherto, and 
now have some individual significance. The probable error is about 
0.02. 
GENERAL DISCUSSION 

The first point of interest is to compare the observed central 
intensities with the values derived from the standard theories of 
Fraunhofer lines. The comparison may be made by plotting r, 
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against W; but if W is in wave-length units, the characteristics of 
the different spectral regions are demonstrated more clearly by 
plotting 7, against W/\. The experimental results are shown in 
Figure 1. The three regions near 6300, 6g00, and 7700 A represent 
the data in Tables 1, 2, and 3; while the six lines near 5550 A are 
taken from the observations of Shane,° whose lines are not quite so 
deep as the author’s. The 4500 A curve is from my earlier work’ and 
may not be reliable. The scatter of the points in this diagram is 
real and represents the individual peculiarities of the Fraunhofer 
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Fic. 1.—Correlation of central intensity with W/A 
lines. We are interested at present, however, only in the shape of 
the mean curves for the various regions, and in drawing these curves 
all lines that are obviously abnormal have been ignored. It is seen 
that the curves differ to some extent for the different regions but 

that they coincide for the weakest lines. 
In order to obtain the corresponding computed relation between 
r, and W/X, it is necessary to know the Doppler broadening, the 
damping factor, and the F-curve (i.e., the curve relating the intensity 
ry at any wave-length within a line to the logarithm of the line-scat- 
tering coefficient at that wave-length). Choosing thermal motion 
corresponding to T = 6000° K and molecular weight 50 for the 
Doppler effect, a damping coefficient of w = 27 X damping breadth 
= 4.2 X 104/X sec (from curve of growth), and the F-curves 
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shown in Figure 2 representing (a) pure scattering,’® (b) Edding- 
ton’s formula™ with € = 0.0, (c) same with e = 0.1, (d) a non- 
theoretical relation, we obtain the A group of curves in Figure 3. 
They are independent of wave-length, and we note that steeper 
F-curves form deeper or sharper lines. Now comparing the computed 
curves with the observed, we find at once the fainter observed lines 
are too wide by a factor of about 1.5. No reasonable F-curve will 
improve the agreement, the damping factor for faint lines has prac- 
tically no influence, and the only way to bring about agreement is 
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Fic. 2.—F-curves: (a) pure scattering, (b) Eddington’s formula with e = 0.0, 
(c) same with e = o.1, (d) a non-theoretical relation. 


to adopt a larger Doppler broadening. This may be done by assum- 
ing that turbulent movements exist in the sun’s atmosphere, the 
velocities of which are comparable with the thermal velocities. The 
consequences of such turbulence on the structure of Fraunhofer 
lines will be considered in some detail. 

In order to study the effect of turbulence on the computed con- 
tours, it is necessary to decide whether the turbulence is on a large 
scale or a small scale in comparison with the depth of the sun’s re- 
versing layer. 

For large-scale turbulence the lines emerging close to the sun’s 
surface have a breadth determined by thermal motion alone, and 


10 A. Schuster, Ap. J., 21, 1, 1905. 1 AS. Eddington, M.N., 89, 620, 1929. 
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the increased broadening that is actually observed in a terrestrial 
spectrograph is due to haphazard displacements of these lines. The 
B group of curves in Figure 3 shows the result of imposing large- 
scale turbulence with a Maxwellian velocity distribution and a 
most probable velocity of 1.8 km/sec. These curves agree fairly 
well with the observed curves and at least give the impression that 
close agreement could be obtained by a suitable choice of the F- 
curve. Such turbulence does not affect the equivalent breadth or 
the curve of growth, and hence the agreement with theory found for 
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Fic. 3.—Calculated relation between r, and W/d for F-curves of Fig. 2, with (A) no 
turbulence and (B) large-scale turbulence of 1.8 km/sec. 


the observed curve of growth”” is still to be expected. One might 
expect, however, that a spectrogram of a well-focused solar image, 
taken with good observing conditions, would show measurable 
velocity shifts. Evershed’ finds that shifts of the order of 0.3 km/sec 
are observable between center and limb but are rarely to be found 
right at the center of the disk. Even the velocity of 0.3 km/sec is too 
small to explain the line broadening, and hence it is at least doubt- 
ful whether the turbulent movements can involve complete granules. 

If, on the other hand, it be assumed that the turbulence is on a 
small scale, the whole contour of the absorption coefficient for each 
line is broadened. Computing the relation between r, and W/) for 


2D. H. Menzel, Ap. J., 84, 462, 1936. 13 M.N., 94, 96, 1934. 
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this case and using a turbulent velocity of 1.5 km/sec, we obtain the 
results shown in Figure 4, which are rather similar to the B-curves 
in Figure 3. The experimental results are in fair agreement with 
both series of curves and do not permit a decision as to which is 
preferable. When the curve of growth is computed with the whole 
(turbulent + thermal) motion, it falls higher on the log W/X axis 
and thus destroys the former agreement between theory and ex- 
periment.” It will be remembered, however, that these com- 
putations were made on the basis of F-curves representing either 
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scale turbulence of 1.5 km/sec. 


Pannekoek’s computations (case B) or pure scattering. If, instead, 
we use the F-curve representing Eddington’s formula, which agrees 
rather better with observations,” "4 the curve of growth is brought 
lower on the W/) axis, and again agrees with observation. Figure 5 
shows that the computed curve of growth based on Eddington’s 
formula (€ = o) anda damping factor w = 5.0 X 104/\ agrees with 
my empirical curve. Thus the hypothesis of small-scale turbulence 
can fit the observations at least as well as large-scale turbulence, 
and at present it is scarcely possible to decide which alternative 
gives the better representation of the actual conditions. 

The breadth of Fraunhofer lines is not the only evidence of 


14M. Minnaert and J. Houtgast, Zs. f. Ap., 12, 81, 1936. 
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turbulence in a stellar atmosphere. Struve and Elvey’’ have found 
it necessary to adopt very large turbulent velocities to explain the 
curve of growth of some early-type stars. In the chromospheric spec- 
trum also, Unséld™ has deduced turbulent velocities as high as 
15 km/sec to explain the structure of Cat lines. Even the appear- 
ance of the sun’s granulated surface would lead one to expect a 
certain amount of turbulence in the reversing layer. The granulated 
areas themselves move with mean velocities of 8 km/sec," and gas 
movements which are an appreciable fraction of this amount must 
be involved. 

















0 2 3 4 5 
log number of absorbing atoms 


Fic. 5.—Curve of growth computed for a turbulent atmosphere and Eddington’s 
formula. The points represent the observed curve. 


It is well to consider briefly what other phenomena might pro- 
duce the observed broadening of Fraunhofer lines. Hyperfine 
structure will widen a few lines and should have much the same 
effect on the curve of growth as small-scale turbulence, but would 
not affect all lines with the regularity shown in Figure 1. The exist- 
ence of noncoherent radiation processes might result in light being 
absorbed in the wings of a line and re-emitted in the center with a 
consequent broadening effect. No adjustment of the F-curve would 
account for this phenomenon, and the simple method employed here 
for computing Fraunhofer contours breaks down. It will be for 


I Ap. J., 79, 409, 1934. 
© Mt. W. Contr., No. 377; Ap. J., 69, 209, 1929. 
17S. Chevalier, Ann. Obs. Z6-sé, 8, 1912. 
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atomic theory to decide whether such a process could produce the 
effect described in this paper; and until such a theory is available, 
the breadth of lines cannot be regarded as a proof of turbulence. 

In order to explain the variation with wave-length shown in 
Figure 1, it is necessary to assume that the slope of the F-curve de- 
creases as one moves toward the red. Indications of such a change 
have been given before,’ the cause of which is probably connected 
with the fluorescent transfer of radiation from the blue to the red 
by cyclic transitions."** The systematic decrease of equivalent 
breadth with wave-length is probably due to the same cause, and 
also the appearance of emission lines in the red of some stellar 
spectra. 

INDIVIDUAL LINES 

It is not difficult to account for the variety in the sharpness of 

individual Fraunhofer lines. For faint lines the atomic weight is the 


TABLE 4 


INFLUENCE OF ATOMIC WEIGHT ON RELATIVE SHARPNESS OF 
FAINT LINES 





ELEMENT 
a oe | 
O Na, Mg K, Ca Ti | Fe Ba 
‘ i | | 7 . aa I ihe 
’ i ~ | 

Atomic weight....... | 10 | 24 40 45 56 137 
Relative sharpness, d | o:7@-| -@.8% | otas | @:¢6e 1.02 1.18 

| 





most important factor, and its influence on the relative sharpness, d, 
may be computed from the Doppler effect. Adopting a reversing 
layer at 6000° K and 1.6 km/sec turbulence, we obtain the sharp- 
ness for faint lines given in Table 4. As the lines become stronger, 
the influence of atomic weight is decreased until for W/X > 40 X 
10~°, d approaches 1.00 for nearly all atoms. Several of the very 
broad lines (i.e., small d) have a low atomic weight. For the great 
majority of lines, however, the atomic weight lies between 40 and 
58, with, consequently, very little effect on relative sharpness; and 
most of the variations are to be attributed to changes in the F- 


8 R. v.d. R. Woolley, M.N., 94, 631, 1934. 
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curve. It will be noticed that there is more variety among the infra- 
red lines than among the visible lines, which is just what would be 
expected if the F-curves have been influenced by fluorescent radia- 
tion transfer. If our explanation is correct, some lines must have 
very flat F-curves. Generally speaking, the damping factor influ- 
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Fic. 6.—Correlation of relative sharpness (ordinate) with excitation potential 
(abscissa). 


ences strong lines only; and for these the shape of the F-curve has 
so much greater influence that it is not possible to determine a damp- 
ing factor from a single observation of central intensity. There may 
be a few cases where the damping factor is so large that the effect on 
faint lines is observable. 

It has already been shown’ that the sharpness depends statistical- 
ly on the excitation potential, and the foregoing considerations 
show that this fact is probably connected with the dependence of 
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the F-curve on the excitation potential. When d is plotted against 
E.P., as in Figure 6, there is a larger scatter of points than can be 
attributed to varying atomic weight. Evidently the sharpness and 
F-curve shape are not governed simply by the excitation potential 
alone, but depend in a more complicated manner on the term struc- 
ture. 

Special attention may be drawn to the structure of the following 
lines: 

a) The K line 7698 A (E.P. 0.0), although situated well within 
the infrared, is sharper than most lines in the visual region. Its un- 
usual sharpness is undoubtedly associated with the fact that it is a 
resonance line; and it is tempting to conclude that all solar reso- 
nance lines have abnormally low central intensities. The observa- 
tions so far published on the Na D lines® and the Ca g line’? support 
this suggestion, but not the Ca* K line. 

b) The lines of the O triplet at 7771 A (E.P. 9.1) have an abnor- 
mally great breadth, due in part to the low atomic weight. Table 4 
shows, however, that 0.72 is the lowest value of d that can be at- 
tributed to atomic weight alone, and the extra broadening required 
to give d a value as low as 0.62 must be associated with very high 
excitation potential and its consequent flat F-curve. The point that 
requires comment is that, if the F-curve is to affect lines as faint as 
these, the change of slope must set in for very weak absorption (as 
for curve d in Fig. 2). The introduction of a finite ¢ into Eddington’s 
formula does not produce deviations in the F-curve which broaden 
such faint lines, and hence does not give a representation of the true 
conditions. 

c) The broad V line at 6216 A (E.P. 0.3) is typical of the ab- 
normality of V lines in the solar spectrum. 

d) The sharpness of Bat 6141 A (E.P. 0.7) is due both to high 
atomic weight and to low excitation potential. The contour is 
possibly affected by a blended Fe line. 

e) All Mg and Si lines” in the lists are greatly broadened by virtue 
of both their low atomic weight and their high excitation potentials. 


19 Redman, M.N., 95, 742, 1935. 
20 T am indebted to Miss C. E. Moore for the identification of the silicon lines. 
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Moreover, Unséld* has shown that the Mg lines have an abnormally 
large damping factor caused by Stark effect. 

f) The breadth of the Na doublet at 6154 A (E.P. 2.1) is probably 
due entirely to low atomic weight. 


Perhaps the main conclusion to be drawn from this work is that 
central intensities can be explained on the simple F-curve method 
of computing line contours, provided a turbulence of 1.6 km/sec 
be assumed to exist in the sun’s reversing layer. 

The author wishes to record thanks to Dr. W. S. Adams, director 
of the Mount Wilson Observatory, for the facilities placed at his 
disposal; to Dr. S. B. Nicholson for his interest and help with the 
apparatus; and to Dr. G. F. W. Mulders for the use of his step- 
wedges. The work was undertaken with the assistance of a Hackett 
Research Studentship granted by the University of Western Aus- 
tralia. 
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NOTES ON SPECTRA OF CLASS Be 
DEAN B. MCLAUGHLIN 


ABSTRACT 

Observations of the relative intensities of components of double emission lines (V/R) 
and of the emission intensities relative to the continuous spectrum (£/C) are given for 
a number of the brighter Be stars, bringing a previous survey (1932) up to date. 

During the twenty-five years of the Ann Arbor program, three stars previously con- 
stant in V/R have begun to vary, and three others, formerly conspicuously variable, 
have become constant. It is suggested that these represent different stages of the his- 
tory of any V/R variable, and that periods of variation alternate with periods of con- 
stancy in most, if not all, Be spectra. The period of recurrence of variability is prob- 
ably a century or more. 

Several years ago the writer’ published a summary of the Ann 
Arbor observations of the spectra of the brighter Be stars. The 
present paper brings the history of a number of these stars up to 
date. During the interval of five years, unexpected interesting and 
suggestive developments have occurred in the spectra of some of 
these objects. Their implications will be discussed at the end of this 
paper. 

In the notes which follow, E/C refers to the estimated intensity 
of an emission line relative to the adjacent continuous spectrum 
outside the wings of the broad underlying absorption. V/R indi- 
cates the ratio of intensity of the violet to the red? component of the 
double emission lines. 


y CASSIOPEIAE (o"50™7, +60°11’), 2™2, Boe 


Extensive studies by the late W. J. S. Lockyer’ and by Clemin- 
shaw‘ demonstrated that the spectrum was almost invariable from 
1911 to 1928. Small variations then set in, and from 1932 to the 
present time the ratio V/R has varied conspicuously. A marked in- 
crease of intensity of the central absorption lines of hydrogen and 

* Pub. Obs. U. Michigan, 4, 175, 1932. 

2“‘Violet” and “red” are used as convenient designations of the components of 
smaller and greater wave-length, respectively. 

3M.N., 93, 362, 619, 1933; 95, 520, 1935; 96, 2, 1935. 


4 Ap. J., 83, 495, 1936. 
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helium was pointed out by the writer.’ These reached maximum 
strength in May, 1936, almost simultaneously with equality of the 
emission components after the red component had been the stronger. 
Between June and August, 1936, these central absorption lines 
faded again. Since then the violet components of the hydrogen, 
Fe 11, and Si 11 emissions have become relatively stronger; but those 
of helium have remained approximately equal to the red compo- 












































y R IQ15.0 1920.0 1925.0 1930.0 1935.0 
<i i ick ii i i is ae ik i i ns da aa a a a aa a (eek a 

2.0Fr y CASSIOPEIAE a 
1.0} —_@_@ @_« _*_°_°__»s_s__ Fé 
0.5- 4 
2.0- 25 ORrionis ra \ 4 
I.O}—_—_________—_ ee L - ae 
0.5+ a aa 4 
2.0F 6 Monocerotis 4 
1.0 ee ° ° — 
0.5L 4 
2.0} ™ Aaquaru 

a | eee . ‘a 
0.5/4 
2.0} aie 

5 —————— 
bie 8” Lacertat ic we 

ee ee a a a ee a a a ce | 








Fic. 1.—Variation-curves of V/R for Be stars. The points depend on means of H8 
and Hy, except for B' Monocerotis, for which the points are based on HB only. Open 
circles depend on only one or two plates, dots on three to nine, and ringed dots on ten 
or more plates. 


nents. The emission structures of all lines have also become narrow- 
er; but those of hydrogen have shown the effect most conspicuously, 
so that Hy is resolvable with difficulty. The spectrum may be about 
to repeat its anomalous behavior of 1934, when the emission lines 
appeared single on the Ann Arbor plates. The curve of variation of 
V/R is plotted in Figure 1. 


W PERSE! (3"29™4, +47°51’), 4%3, Bse 
One spectrogram, taken August 11, 1932, shows the red compo- 
nents of HB and Hy so much stronger as to be practically conclusive 


5 [bid., 84, 235, 1936. 
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evidence of a temporary variation. Several other plates taken since 
then show V/R nearly 1.0. This star is shifted from the constant to 
the suspected class. 


23 TAURI= MEROPE (354074, +23°39’), 4%3, Bse 
This star, previously variable in both E/C and V/R, has shown 
no marked activity since 1931. 


X PERSE! (3'49™1, +3045’), 6%2-6"9, Boe 


The following observations have been obtained since the earlier 








survey: 
— V/R No. of 
Plates 
| 
‘rf: eae 1.4 4 
Oe, o.oo ae 2 3 
OCT, are 0.8 4 
NG cte c <wiste one oa 2 2 
1930.9... ae 0.7 4 








These data suggest a period a little longer than two years, but it does 
not satisfy the earlier observations. The variation is probably ir- 
regular. 
11 CAMELOPARDALIS (4"57™4, +58°50’), 5™3, Bze 
This star is one of the best examples of variable E/C. Recent data 
for several lines are given in the following table: 





E/( 
DATE ee ee 
| PLATES 
Hy Hé He Ht 

BOBO vigse cele <i eres 0.75 0.4 oe Ga cae a ier 3 
a Se ee vr «2 8 5 0.3 3 
Bi Ping c e oe. A Es 8 7 .4 3 
RO AS 66 ood eer ees 1.6 9 6 4 4 
FOSO6 Ginc ain, ce 2, eGR ¥.6 0.95 o.7 0.5 4 

















The Ann Arbor data for Hy are plotted in Figure 2. Minima of in- 
tensity were reached about 1920 and 1933. The emission has not 
yet recovered the maximum strength exhibited about in 1929. On 
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a plate taken December 14, 1936, the lines Hn, HO, and Hi have 
distinct but faint central absorption, which has never been detected 
before on Ann Arbor plates, though it has been carefully sought. 
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Fic. 2.—Variation-curves of E/C for Be stars. The points depend on 8, except 
for 31 Pegasi, for which curves of H8 and Hy both are given. Open circles represent 
single plates, dots depend on two to nine plates, and ringed dots on ten or more plates. 


105 TAURI (5'02™o, +21°34’), 6%o, B3e 


This star has fulfilled the expectation that it would prove to have 
variable V/R. The complete Ann Arbor data are as follows: 








No. of 

Date V/R Plates 
1630.0. ..... 0.6 2 
1931.8 0.6 4 
1932.8 0.65 I 
oe ee 0.8 2 
1934.8 72 I 
1935.8 1.6 I 
1936.9 Ps I 











The cycle of variation is evidently several years in length. 
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25 ORIONIS (5*19™6, +1°45’), 4¥7, B3e 
A detailed discussion of the variations of this spectrum is given 
by Miss Dodson.® Observations subsequent to her study are given 
in the following table: 











| V/R 
No. OF 
sina PLATES 
HB Hy 

1933.8 1.05 1.0 3 
1044.68... 2 5.52. 0.95 0.9 3 
1934.8 0.75 0.8 4 
BOBS Dik cco a ties wieeies 0.9 1.0: 3 
1935.9 0.8 1.1 4 
$0905:26..66%: 1.0 1.0 4 
ROBO) oo ide os Hares 0.9 1.0 3 














The complete curve of V/R is plotted in Figure 1. The points repre- 
sent means of HB and Hy. Conspicuous changes of V/R continued 
from 1915 to 1932, in cycles varying in length from about 1800 to 
1000 days. Since 1932, however, the star has been quiescent, with 
equal emission components. 


120 TAURI (5'27™6, +18°29'), 5“5, B3e 


The complete Ann Arbor observations are as follows: 




















V/R 
DATE hie 
: PLATES 
HB Hy 

oe. ae 2 I.2 2 
TOGE Se 552 r.2 ee 3 
Oo OC Ae 0.8 0.8 2 
BOBS Geri  ie ces 0.8 1.0 2 
cc? Oe: : Er. 1.1 I 
TOPS Oeics te seke 0.8 0.9 2 
WOAGQ, Wece sce chs ae 2. || I 





This spectrum is probably variable in V/R, with small range. A 
period of about three years is a possibility. 


6 Tbid., p. 180. 
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HD 41335 (5"59"4, —6°42’), 51, Bze 
In the previous note on this spectrum’ the words “red” and 
“violet” should be interchanged. On plates taken in 1926, 1930, and 
every year thereafter until the present, the red components are 
much the stronger, though there has been a distinct change toward 
equality with probable secondary variations superimposed. 


HD 44458 (6"16™8, —11°44’), 5"5, B2e 
The data on this spectrum are fragmentary and indicate simply 
that V/R variation of moderate range has continued. The length 
of the cycle is still uncertain and is perhaps irregular. 


B' MonocEROTIS (624™0, — 6°58’), 4"7, B2e 
Recent Ann Arbor observations are given in the following table: 











V/R 
No. oF 
DATE 
PLATES 
| HB Hy 

200 0): Se ere 2.0 Ls 3 
1934.2 1.6 r.3 3 
1934.5 Py 5 Ae 3 
BO ck as oss Ris t.% 2 
1935.8 I.I 0.8 2 
1936.3 0.8 0.75 2 
1930.9 0.5 0.6 2 














In the study of this spectrum by Miss Hawes* it was shown that 
V/R remained constant at 1.0 from 1905 to perhaps 1923. Between 
that time and 1928.8, V/R decreased to 0.5, and during the next 
two years the inequality was reversed. The foregoing table shows 
that the second reversal in 1935 occurred slowly, and that Hy clearly 
preceded HB, as Miss Dodson has shown in the case of 25 Orionis. 
A curve of V/R for 8' Monocerotis is given in Figure 1, for the inter- 


val 1910-1937. 
x DRACONIS (1229™2, +70°20’), 3“9, Bse 
Periods of thirty and twenty-three years have been suggested by 
Struve® and Jessup,” respectively, for the variation of E/C. Emis- 


7 Pub. Obs. U. Michigan, 4, 182, 1932. 9 Pop. Asir., 33, 596, 1925. 
8 Pub. Vassar Obs., 4, 29, 1934. 10 Ad. J., 76, 75, 1932. 
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sion has completely disappeared at H®, and only a faint trace of 
bright Ha is seen on wide absorption on a plate taken April 16, 
1936. 
20 VULPECULAE (20"07™8, +26°11’), 5“9, B8e 
From about 1921 to 1932 the H® emission increased through a 
small range. In 1936 its intensity was again about equal to that 
observed in 1921. 


25 VULPECULAE (20"17™7, +24°07’), 5™4, B8e 

The HB emission appears to have faded slightly between 1931 

and 1936. 
 CYGNI (20543™5, +36°07’), 4%5, Bse 

Weak single emission was present at H6 in 1926. In 1931 there 
was no emission, and plates taken each year after that showed HB 
as a normal absorption line with a rather narrow central core. Four 
plates taken in 1936 show faint emission once more, and the absorp- 
tion appears wide and ill defined, as in 1926. 


HD 199356 (20"51™6, +39°51'), 7Mo, Boe 
There are only two spectrograms of this star, but they definitely 


establish the variation of V/R. 


Date V/R 
$088 Ae, Gi... ens ..0.6 
0G ol) an re 2 


ft CyGni (20"56™4, +47°08’), 49, Boe 
This interesting star lost its emission lines in 1912, regained them 
faintly in 1925 and strongly in 1927, and showed conspicuous varia- 
tions of V/R in 1928-1930. The HB emission faded again until it 
was just perceptible above the continuous spectrum in 1932. From 
that time to 1936 it remained quite inactive, without noticeable 
loss or gain of intensity. 


uv CYGNI (215138, 34°20’), 4"4, B3e 


This star shows definite V/R variation of small range. During 
1929 the emission components were approximately equal. From 
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1930 to 1933 the red components were the stronger, and since 1934 
the violet components have been slightly more conspicuous. If there 
is a period, it is of the order of eight or ten years. 


6 CEPHEI (21517™3, +64°27’), 5"2, B3e 
From 1932 to 1936, E/C at Hy changed from 1.0 to 0.75. A 
previous minimum of intensity, 0.7, occurred in 1921-1924. In 1936 
the red component of Hy was slightly the stronger, though in previ- 
ous years the components were practically equal. 


€ CAPRICORNI (21531™5, —19°54’), 4"7, Bsep 
This spectrum is characterized by well-defined hydrogen emis- 
sion with strong sharp central absorption, and by variable enhanced 
metallic absorption lines. The latter emerged to visibility in 1929, 
faded in 1931, remained absent through 1935, and reappeared previ- 
ous to September, 1936. 


o Aquaril (21558™1, — 2°38’), 4"7, Bse 
Twenty-seven plates taken between 1932 and 1936 show no indi- 
cation of variation of either E/C or V/R. 
25 PEGASI (2203™1, +21°13’), 5"7, Bse 


The following observations have been obtained: 








Date E/C (HB) No. of Plates 
a | 0.5 2 
1043-0). ...« 4 I 
1934.8 : 5S= I 
re re — I 
FOSO 0.6 6% 4% 0.25 I 











The gradual fading of the emission previously noticed, between 1917 
and 1931, has evidently continued. The curve of variation is shown 
in Figure 2. 
31 PEGASI (22"16™6, +11°42’), 4M9, B3e 
Recent Ann Arbor observations are summarized in the following 
table: 
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E/t | 
DaTE eee ata | eat 
| PLATES 
Hg Hy 
1932.7 25 | 1.3 | 6 
1933.6 3-0 | I.4 5 
1934.6 2.5 is3 | 4 
1935.6 1.0 0.4 4 
1936.6 0.4 | 6 


Between 1933 and 1936 the emission faded from the brightest to 
the faintest ever recorded in this spectrum. A curve of variation 
is shown for HB and Hy in Figure 2. A previous minimum of emis- 
sion occurred between 1912 and 1919. 


mt AQUARII (2220™2, +0°52’), 4™6, Bie 
Ann Arbor observations of V/R between 1911 and 1930 are tabu- 
lated in a previous paper.’ During that interval the value of £/C 
at HB was 2.3, with possible small fluctuations. Later observations 
of V/Rand E/C are given in the following table: 




















“| 
DATE F B/C signed 
(HB) PLATES 
HB Hy | 
1931.4 2.8 "5 2.0 I 
1931.55 1.05 1.0 5 4 
1931.05 0.44 0.62 1.4 5 
1931.80 0.36 °.40 1.6 5 
1932.45 1.5 1.5 1.3 5 
1932.0 2.9 1.0 5.2 0 
1932.8 7 1.8 bd 5 
1933.4 I.4 1.0 t.2 4 
1933.0. r.0 1.0 3.3 4 
1933.8 0.8 0.9 re 2 
IOG4.4.. :: 0.0 0.0 I.4 2 
TOS4: O's. 5: 0.4 0.4 1.6 4 
1934-75 0.3 0.4 2.0 4 
TOG: Ss. « I: ES I.1 3 
1935.0. B.S I.4 ee 5 
BAGO yo avsnon 1.0 °.9 1.0 4 
Uo se Ry eee 0.75: . 1.0 3 
Cr iy eee ae : a3 0.7? 7 














1 Pub. Obs. U. Michigan, 4, 44, 1931. 
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The recent variations of V/R have been complicated by the fading 
of the emission, which rendered the numerical values unreliable 
although the direction of the inequality is definite enough. The 
data are in satisfactory accord with those given by Lockyer.” 

The curve of V/R from 1911 to 1936 is shown in Figure 1. Un- 
less the plates earlier than 1924 were taken, by remote chance, at 
times of equality of emission, the ratio V/R was constant near 1.0 
from 1911 to 1923, and the large range of V/R is a comparatively 
recent development. This case is thus similar to those of y Cassio- 
peiae and 6" Monocerotis. 

Another point of interest is shown by the tabular values of V/R 
in 1932-1933. Hy passed through equality of components a few 
months before H@. This is similar to the behavior of these lines in 
25 Orionis and 8' Monocerotis. 


8 LACERTAE (NorTH) (22"31™4, +39°07’), 5™8, B3e 
Observations obtained since the previous survey are given in the 
following table: 

















V/R 
DaTE No. oF 
PLATES 
Hp Hy 
£032 07.55 1.3 OR 4 
1933.0.. 2 2 4 
1934.7 2.8 2.2 3 
1935-7 | 0.9 0.8 3 
1930.7 0.5 0.7 6 
| 








The curve of variation of V/R is shown in Figure 1. A period of 
about eleven years appears to be a possibility. 


HD 217050 (22"52™7, +48°09’), 5“2, B3ep 
This star is a good example of a variable composite spectrum. 
Emission lines were absent in 1921. Since 1928 they have been 
strong, with sharp strong central absorption. With them an en- 
hanced metallic absorption spectrum appeared, and it has persisted 
unaltered through 1936. 


2 M.N., 96, 480, 1936. 








































NOTES ON SPECTRA OF CLASS Be IgI 


8B Pisctum (22'58™8, +3°17’), 4"6, Bse 
This star shows conspicuous £/C variation. Observations since 
the previous note are given in the following table: 








E/¢ 
DATE ———_________— sind No. OF 
| PLATES 
Hp Hy Hs 
OC ree 0.7 0.3 4 
1933-5 1.1 5 0.25 | 4 
Ma eee I.4 7 3 4 
1034-7---++--- 2.2 .85 5 6 
a > ae 2.1 8 5 4 
ies Eh eer ere 2.0 0.8 | O.: 4 
| 














The curve of variation is shown in Figure 2. In view of the relative- 
ly short duration of the recent minimum, observations about 1915 
would be of considerable interest. 
DISCUSSION 

Several stars which have shown no new developments have been 
omitted from the preceding notes. The recent observations have re- 
sulted in a change of status of eleven stars. Some merely represent 
confirmation of suspected variation; others represent radical changes 
of behavior. The eleven stars are listed in the following table, with 
their classification as of 1932 and the revised classification of 1937. 

















V/R E/( 
STAR —s — 

1932 1937 1932 1937 
2 Sk | ae eae re const? var!! const! const! 
OH, .. rnreee Ay const var? const! const! 
WOR BOM isc oes var! var! Zecoae const 
ra |S i ee ae eee var!! var!! const! var? 
ot ee var | ne See ae const 
ie ee var!! var!! const 
jc 07 ee Se var!! var!! wre const 
re TR see const? const? var? var 
a0 Wille e ss 65% const const const var 
OLR so asic cas const? var? var? var 
MN a a cuees var!! var!! const! var! 

















In addition to these, the star HD 199356, not included in the previ- 
ous survey, has been found variable in V/R. 
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In the previous survey of forty-three Be stars (omitting the pe- 
culiar cases 8 Lyrae and v Sagittarii), it was concluded that twenty 
showed variable V/R, eight were more or less strongly suspected of 
variation, six were certainly constant, and nine were uncertain be- 
cause of insufficient data or single emission. Now, twenty-one are 
recognized as variable (y Cassiopeiae added), eight are suspected 
(y Cassiopeiae confirmed, wy Persei added), five are constant, and 
the same nine are still unknown. 

In 1932 the statistics of E/C were: twelve variable, three suspect- 
ed, twenty constant, and eight unknown. These figures are now 
modified to: sixteen variable, two suspected, twenty-one constant, 
and four unknown. 

The general conclusion that variability is the rule and constancy 
the exception among Be stars is slightly strengthened. The conclu- 
sion that no stars with large range of V/R show change of E/C is 
now disputed by the fading of emission in 7 Aquarii. It now appears 
more probable that it is only a question of time before all Be stars 
will show variable E/C. The irregularity of lengths of cycles serves 
to emphasize the earlier conclusion that binary motion as an explana- 
tion of the changes must be ruled out in most cases. 

The present survey and work done, since the previous survey, 
on the earlier history of some of the stars has revealed cases of 
cessation of variation of V/R, as well as some cases of incidence of 
such variation in stars previously constant. Thus, on the one hand, 
we have ¢ Tauri,’ b? Cygni,"4 and 25 Orionis, which formerly showed 
large changes of V/R but are now constant. It may be mere coin- 
cidence that the last large inequality of components shown by all 
three of these stars was in the same direction, namely, V/R greater 
than 1.0. 

On the other hand, three excellent examples of incidence of V/R 
variation are now recognized. These are: y Cassiopeiae, 6* Mon- 
ocerotis, and w Aquarii. These do not present the same unanimity 
with regard to the direction of their first departure from V/R equal 
to 1.0 as did the previous three in their last large inequality. 

It is possibly a coincidence that the numbers of cases of incidence 

13 Losh, Pub. Obs. U. Michigan, 4, 3, 1931; McLaughlin, zbid., 4, 181, 1932. 


14 Losh, tbid., 4, 200, 1932. 
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and cessation are equal, but it has at least the suggestion of statisti- 
cal equilibrium of the class. To the writer it is inconceivable that 
these observed changes represent permanent changes in the status of 
these stars. Rather, it seems very plausible that a given star may 
have alternate periods of activity and quiescence. In y Cassiopeiae 
we have been privileged to observe the transition in one direction; 
in 25 Orionis the reverse change has been observed in detail. In the 
course of time we may hope to see the resumption of variation in the 
spectrum of 25 Orionis and its cessation in y Cassiopeiae. When we 
consider the microscopic size of our time interval of observation com- 
pared with the life of a star, it would be illogical in the extreme to 
postulate a unidirectional change in any individual case when we 
observe equal numbers of cases in which the change has occurred in 
two opposite directions. 

The present survey, added to the earlier work, shows that the 
periods of variation of £/C are, for the most part, many years in 
length, with ten years as a probable lower limit. Periods of two to 
twelve years would include almost all of the V/R variables, when 
we consider only the cycles of V/R. The periods of recurrence of 
variation or constancy of V/R are perhaps of the order of a century 
or more in length. 

The fact that any given Be spectrum has shown no variation over 
an interval of twenty years or so is no guaranty that it will remain 
forever an uninteresting object. The case of y Cassiopeiae is suffi- 
ciently recent to serve as a reminder that every Be star is at least 
potentially variable, and hence deserves reobservation at intervals 
sufficiently short to give warning of approaching activity. Only 
when we have a larger number of observations of incidence and cessa- 
tion of variation will it be possible to generalize confidently with re- 
gard to the manner of such incidence or cessation, and only then will 
it be possible to attack the problem of physical causes with any con- 
siderable hope of success. 

THE OBSERVATORY 
UNIVERSITY OF MICHIGAN 
February 27, 1937 








ON THE INTERPRETATION OF THE SURFACE 
BRIGHTNESSES OF DIFFUSE 
GALACTIC NEBULAE 


OTTO STRUVE 


ABSTRACT 

Observational results indicate that the surface brightnesses of many dark nebulae 
are not uniform but exhibit dark cores and more luminous edges. This effect cannot 
be due to uneven thickness or transparency of the dark nebula. The Lommel-Seeliger 
theory of diffuse reflection is extended to the case of an opaque spherical nebula consist- 
ing of particles which are characterized by different spherical albedos and phase func- 
tions. The surface brightness of the nebula is computed for different distances from the 
center, under the assumption that the illumination comes from the uniform brightness 
of the starlit sky. The results give values for the difference in stellar magnitudes be- 
tween edge and center, from 0.03 mag. to 0.17 mag. The observations favor slightly 
a large albedo and a phase function which throws most of the light in the direction of 
the incident beam. The Lambert phase function can only be reconciled with the ob- 
servations if the albedo is assumed to be unreasonably large. 

The Lommel-Seeliger theory is next applied to the peculiar phenomenon of narrow 
bright rims frequently observed in emission nebulae but never found in nebulae having 
continuous spectra. A list of these objects is given, the most conspicuous of which is 
the rim on the west side of the Horse-head nebula, Barnard 33. If in the derivation 
of the Lommel-Seeliger formula the absorption coefficient of the incident radiation is f; 
and the absorption coefficient of the visible diffuse nebular light is #2, then for a nebula 
in = the product of the radius by &; is very large, the surface brightness of the 
nebula is 


cos 7 
cos € + (k2/kr) cosi’ 





I = const 


where 7 is the angle of incidence and e is the angle between the line of sight and the 
perpendicular to the nebula. If &2/f; is negligible, J is inversely proportional to cos e, 
so that at the edge a narrow bright rim will be formed. 

1. The principal purpose of a series of investigations now being 
carried on at the Yerkes Observatory is to obtain data concerning 
the nature of the particles which give rise to the scattered light of the 
nebulae. Henyey’s recent paper’ on this subject has elucidated the 
theory of nebular scattering and has shown how, with reasonable 
assumptions concerning the form of the phase function, the geometri- 
cal conditions of the problem and the albedos of the particles may 
be derived. From the observations of the colors of a few of the nebu- 
lae and from the small amount of radial polarization observed in 
two or three of them, we conclude that pure Rayleigh scattering de- 


* Ap. J., 85, 107, 1937. 
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pending upon A~ can, at most, contribute only a fraction—perhaps 
not more than one-half or one-third—to the total light of these re- 
flection nebulae. It is, however, certain that some nebulae show se- 
lective reddening for transmitted light, so that scattering by small 
particles is not altogether negligible. The observations of the colors 
of the nebulae, of the amount of selective reddening for transmitted 
light, and of polarization really prove only that the majority of par- 
ticles are not small enough to scatter in accordance with the sym- 
metrical phase function g(a) = 1 + cos*a. The particles either may 
be large enough to have phase functions which return most of the 
light toward the source or may be sufficiently small to have phase 
functions which throw most of the light in the direction of the inci- 
dent beam. 

While Henyey and the writer have used the Lambert phase func- 
tion without having encountered serious inconsistencies in the re- 
sults, other phase functions would have served equally well. The 
trouble lies in the bewildering structural formations which are pres- 
ent in practically all diffuse nebulae and which tend to obliterate the 
effects produced by the geometrical conditions and by the form of the 
phase function. It will evidently be a difficult task to determine, 
even roughly, the form of the phase function. Yet, this is required 
if we are to make progress in our study of the physical properties of 
the nebular particles. 

Since there are apparently no simple cases of nebulae where the 
phase function can be derived directly from the surface brightnesses 
at varying distances from the illuminating star, it will be necessary 
to pursue the laborious method of separating and interpreting those 
nebular features which are not merely a result of density fluctua- 
tions. The present investigation makes a first attempt in this direc- 
tion. We shall show that while, even now, it is not possible to derive 
the phase function accurately, a number of interesting features, com- 
mon to many nebulae, can be satisfactorily explained by the theory 
of diffuse scattering. 

2. In a recent paper? it was shown that the dark nebula Barnard 
15 has a dark core and slightly more luminous edges, although in 
all parts of the nebula the faint stars are completely blotted out. 


2 Struve and Elvey, ibid., 83, 172, 1936. 





196 OTTO STRUVE 


The total range is only o.1 mag. between center and edge, and the 
observations are satisfactorily represented by a straight line. Photo- 
graphs of Brs show conclusively that this effect is not produced by 
the increasing thickness of the nebula as we pass from the edge 
toward the center. Even at the adopted edge, marked by abscissa 
1.0 in Figure 1, the absorption of starlight behind the nebula is still 
of the order of several magnitudes. Hence, the background of the 
starlit sky behind the nebula is effectively blocked out. 
The photographs of the Milky Way by Barnard and by Ross re- 
veal a large number of dark nebulae possessing similar dark cores. 
Mag. 
11.47 





11.49 - 


II.51 






Region of transition 
Star. cloud 








foes Edge of opaque nebula 


11.59 ! | l | | 
° 0.2 O04 0.6 #08 1.0 Distance from center 





Fic. 1 


Plate IV shows the extensive dark nebulae in Taurus, reproduced by 
E. L. McCarthy by superposing two original negatives taken by 
Professor Ross. The dark nebulae are very conspicuous, and the 
cores are quite prominent. Plate V shows the same region, repro- 
duced from an original negative taken by Professor Barnard with 
the 6-inch Bruce telescope. The black cores are almost identical in 
the two illustrations. Finally, Plate VI is a reproduction obtained 
from the superposition of two original negatives by Barnard: it il- 
lustrates the dark nebulosities in the region of a 4"10™ 6 + 54°30’. 
Barnard’s dark nebulae B11, B12, and B13 have well-marked black 
cores; but a large starless region, B8, has a surface brightness which 
exceeds that of the cores and resembles that of the background in 
the denser star clouds. This region shows very clearly that there are 
marked differences in the surface brightnesses of dark nebulae. 
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PLATE V 
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BARNARD’S DARK NEBULAE IN TAURUS 
From a photograph taken with the 6-inch Bruce telescope 


(The numbers of Barnard’s dark nebulae are shown on the margins) 
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The color of the luminous edges surrounding the cores of dark 
nebulae is yellow. This was shown on photographs of the dark region 
north of 22 Scorpii, photographed with the Schmidt camera in yel- 
low and in blue light.’ More recent photographs of Barnard 216 
confirm this result. The marking described by Barnard in his Cata- 
logue of 349 Dark Objects in the Sky as a “dark spot in lane” has an 
appreciably smaller surface brightness than the surrounding starless 
region. It is equally well visible on an Eastman Process plate ex- 
posed without a filter and on an Eastman IC plate exposed with a 
red filter (Wratten 23 A). We conclude that the luminous portions 
of dark nebulae resemble in color the background of the night sky, 
and are therefore similar to the integrated light of the stars. 

3. Since reasons have already been given for believing that an 
opaque mass of scattering particles in interstellar space illuminated 
by general starlight will be bright enough to produce a measurable 
amount of surface brightness on top of the luminous background of 
the sky, we shall follow up this hypothesis more closely and investi- 
gate whether it is possible to account for the uneven surface bright- 
ness of the dark nebulae. 

Let us first consider the dark nebula as an opaque spherical body 
which acts as a diffusing screen obeying the law of Lambert. If the 
incident light is designated by L, we find for the surface brightness, 
per steradian, of the nebula 


A ; 
I =— Locos?. 
us 
Here A is the Lambert albedo. If we imagine that Z is uniformly 
spread over one hemisphere so that L, is the surface brightness per 


square degree of the sky, as seen from the nebula, such that 


VBP tide we b, 
40 


we obtain for the surface brightness of the nebula per square degree 
A me ae 
IL =—L, 27 cosisinidi = AL,. 
Tw ° 


3 [bid., 84, 223, 1936. Recent photographs in red light confirm this. 
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L, is not accurately known. The best we can do now is to use the 
value of {Z, dw given by Seares, Van Rhijn, Joyner, and Rich- 
mond‘ for the total amount of starlight: 


285 {, dw = 577 stars of photographic magnitude 1.0. 
Hence, we find for one square degree of the sky a surface brightness 
of 56 stars of magnitude 10. This agrees quite closely with the re- 
sults of Bottlinger.s The nebula receives, of course, more light be- 
cause the units we have adopted refer to starlight which has passed 
through the atmosphere. Since, however, we observe the nebula also 
through the atmosphere, we can disregard the corresponding cor- 
rection. 
If we adopt A = 0.3, we find 


I; = 17 stars of magnitude ro per square degree. 


This amount of light is superposed over the “‘terrestrial’’ background 
of the sky, for which we shall adopt a value of 4.4 pg. mag. per square 
degree, which is equivalent to 174 stars of magnitude 10. This value, 
due to Dufay,°® must be corrected for the light of unresolved stars 
fainter than those which could be eliminated in his observations. 
This amount may be estimated at 20 per cent, leaving the equiva- 
lent of 140 stars per square degree for the “‘terrestrial’’ background. 
Consequently, if Lambert’s law is adopted and if A = 0.3, the nebu- 
la should appear as a surface which is slightly brighter than the 
“terrestrial” sky background, the difference being 


Amr, = 2.5 log 13} = 0.12 mag. 


While this value is quite plausible, Lambert’s law leads to an ex- 

pression for J; which is independent of the angle ¢ between the line 

of sight and the perpendicular to the diffuse surface. Accordingly, 

I, is constant for all €, and there should be no variation in the sur- 

face brightness of a dark nebula. It is not likely that this inconsist- 
4 [bid., 62, 373, 1925. 


SZs.f. Ap., 4s 373, 1932. 
© Etude de la lumiére du fond du ciel nocturne, Paris, 1934. 
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ency with the observations can be removed by dropping our simplify- 
ing assumption that the starlit sky, as seen from the nebula, appears 
uniformly bright and that its 
surface brightness is equal to the p Lilie 
average surface brightness of all 

stars as seen from the earth. 
Hence, we are compelled to 
abandon Lambert’s law. 

4. We shall next use the law 
of diffuse reflection proposed by 
Lommel and by Seeliger. If the nebula is a uniform plane-parallel 
slab of diffuse material having an absorption coefficient k and a 
thickness R, the surface brightness is’ 








FIG. 2 


cos t 
a 


I =* Lea) —— 


oe e— Rk(sec i+sec ©] 
cos i + COS € 


The factor ng(a) measures the fraction of the incident light which 
is reflected by unit volume in the direction of the phase angle a. If 
the nebula is opaque, so that Rk is large, we have 


, cos 7 
I="L —, I 

o(a eT ae (1) 
No assumption has here been made concerning the nature of the 
nebular particles. If, however, we assume that we are dealing with 
fairly large spherical particles of radius p, we have 


k = Nap’, 
a P 
w= Tp’, 
27 


fe(a)dw = 27. 


The albedo y gives the fraction of the total light incident upon each 
particle which is reflected in all directions. Then 


Cos 2 


ee a 
.* aT Lo(a) cos i + cose’ 


7 Schoenberg, Handbuch der Astrophysik, 2, Part I, 35, 1920. 
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paper,® where y = 2p. If we again replace Z by L,dw and integrate 
over the hemisphere which is visible from the nebula, we obtain the 
surface-brightness J, per square degree of the nebula as seen from 
the angle e: 


COs 1 


ae = 
21 Lf oe cos 1+ cos € al (2) 


Since the phase function is not known, we shall make three as- 
sumptions for it and compute the surface brightness of a nebula 
which is uniformly illuminated from all sides by starlight. The in- 
tegration will therefore be carried out for six cases, as shown in the 
accompanying table. 








PHASE FUNCTION 








alls 5) nal 
er(a)=Cos” > | ¢2(a)=0.5 g(a) =sin* = 
3 CG naa ree ee (1) €=0 | (3) e=o (5) e=o 
OS Ce ae (2) e€=90° (4) €=90° (6) e=go° 





Case (1) is identical with the one computed in our former paper: 
Is = 0.2570, . 


Our phase function ¢,(a) is that of Euler, which agrees reasonably 
well with the one by Lambert and throws most of the light in the 
direction opposite that of the incident beam. ¢;(a) throws most of 
the light in the direction of the incident beam; and while its form 
has no physical justification, it resembles the phase functions of 
moderately small particles computed by Blumer from the theory of 
Mie. ¢,(a) is an intermediate case. We shall assume that formula 
(2) holds for all three expressions of g(a), although it was derived 
from the more general expression (1) under the assumption that the 
particles are large enough to permit us to use k = Nap’. However, 
if we compare (1) with (2), we find that their forms are identical 


8 Ap. J., 83, 164, 1936. 


which is identical with Seeliger’s expression adopted in our former 
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and that only y changes its meaning if the particles are so small 
that diffraction effects became serious. 

It is obvious from the symmetrical shape of the three phase func- 
tions and from the identities 


t1¢r(a) + gr(a + 180°)] = 0.5, 
r1¢3(a) + 9;(a + 180°)] = 0.5, 
$2(a) sa 0.5 ’ 
that at the edge (e€ = go°), J, will be the same for all three functions. 
The integrations are carried out as follows: 
Case 1: See Astrophysical Journal, 83, 166, 1936. 
Case 2: Let 8 be the azimuth of the source as seen from the nebu- 
la. Then 


cos a = sinicos#, 


¢i(a) = 0.5 + 0.5 sinicos®B. 





The integral becomes, since ——._——— = 1, 
cos 2 + COS € 


“w/2 2m 
[= cf { ¢:(a) sini didB = Cr. 


Case 4: Similar to Case 2. 
Case 6: Similar to Case 
Case 3: € =o;t=a 
@/2 
sin acosa 
I = 0.5C anf ee de 
° 


1+ cosa 


¥ 


NO 


Substituting 1 + cosa = 2, 
72-1 
I =0.5C anf —— dz = 3 — Inz|{ = 0.307Cr . 
I ~ 


Case 5:€ =o;i=a 
W/2 
P @ . 
sin? - sin acosa 


I=C 2 ih one da. 
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ae } 
Substituting 
4 
° a 
| sin’ > = 0.5 — 0.§ cosa, 
I+ cosa =z, 
( 2 g- 2 g-— 2 ) f 
I= Carlos =+ ds —0.5f ( — ssi Cita { 
\ I > I ad : 
f 


Table 1 gives a summary of the results. 


TABLE 1 





















Phase Euler: 
Function ¢1(a) =cos? 4 
€=0 Is=0.25yLi 
€=9g0 Is=0.57L 
Ratio: 
Edge/Center 2.0 


. Sma articles: 
Constant: mall Particle 


- a 
¢2(a)=0.5 ¢3(a) =sin? : 


Is3=0.15yLy Iss=0.06yL; 


Isg=0.5yLr i6=0.5yLi 
i Rs SA RY (Pha ae Ss 3 : 
3-3 8.3 
























Assuming y = 0.3 and L, 


degree, 


Tax I; 


ror 


= 8.4, fe 


For ¢,(a): y = 0.3; Am 


For ¢.(a): y 0.3; Am 


For ¢;(a): y = 0.3; Am 





If we again adopt for the “‘terrestrial’’ background of the sky the 
equivalent of 140 stars of magnitude 10 per square degree, we ob- 
tain for the difference in magnitude between the edge and the center 
of a dark nebula the following values: 


t4 I,, = 1.0, 
oe Io = 8.4 
; 
45.4 ; 
2.5] = 0. ag. 
5 log ong 0.03 mag 
8. 
2.5 log "ws 0.04 
142.5 
8. 
2.5 log 4°:4 = 0.06 





DIFFUSE GALACTIC NEBULAE 203 


The observed value is 0.1 mag., or twice the result computed for 
y;(a). Increasing the value of y to 0.6, which would correspond to a 
rather efficient scattering power of the nebula, we obtain: 


For ¢.(a): y = 0.6; Am = 0.06 mag. 


For ¢.(a): ¥ 0.6; Am = 0.08 


For 9,(a): y = 0.6; Am = 0.11 
The observations seem to favor the larger value of y and the phase 
function ¢,(a), which throws most of the light in the direction of the 
incident beam. It is possible, however, that y is even larger, in which 
case ¢.(a) or g,(a) would give plausible values for Am. Thus, for the 
unreasonably large value of y = 1.0, we find: 
For g:(a): y = 1.0; Am = 0.09 mag. 


For g,(a): y = 1.0; Am = 0.13 


For g;(a): y = 1.0; Am = 0.17 


Our result must be regarded as tentative because it is based upon 
uncertain values of the ‘terrestrial’ background and of the real 
amount of starlight as seen from the nebula. Also, the observations 
are somewhat uncertain and should be repeated for a greater number 
of objects. Nevertheless, it is not very probable that we were able, 
with the Fabry photometer, to detect the full range of the variation 
in surface brightness. The computations give us the maximum range. 
Hence, our conclusions regarding y and ¢(a) are probably correct. 

5. In the derivation of the Lommel-Seeliger formula for diffuse 
reflection, the effect of the reflections of orders higher than the first 
is neglected. It is of interest to find to what extent the foregoing 
conclusions will be affected, if the reflections of the second order are 
allowed for. This can be done very simply by means of the formulae 
for diffuse reflection derived by Fessenkoff? and by Schoenberg.” 

Fessenkoff adopts for the phase function Rayleigh’s law: 


g(a) = 1 + cosa, 


9 Bull. Soc. Astr. de Russie, May, 1916. 10 Op. cit., p. 37. 
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and derives an expression for the surface brightness, taking account 
of the second-order reflections. Schoenberg uses 


y;(a) = 1 — 2.7 cosa + 3 cCos?a, 


which is the phase function for small drops of water resulting from 
Mie’s theory. ¢,(a) is symmetrical with respect to both the line of 
sight and any straight line passing through the particle at right 
angles to the line of sight. Hence it should give results which are 
comparable to our ¢.(a). On the other hand, ¢,(a) throws most of 
the light in the direction of the incident beam, and is therefore some- 
what similar to our ¢,(a). 

Schoenberg has computed tables for both formulae, giving the 
value of the surface brightness of the diffusing screen at intervals 
of 10° in 7 and in ¢, and at suitable intervals of the azimuth A. The 
integration for the hemisphere can therefore be carried out numeri- 
cally: we first compute from the tables of Schoenberg the average 
surface brightness for each value of 7. This eliminates the azimuth. 
Next we multiply each surface brightness thus computed by the 
area of the corresponding zone in 7. The summation of the results 
for all zones gives the relative surface brightness of the nebula 
for any given e. Setting the surface brightness at ¢ = o equal to 1, 








we have 
TABLE 2 
. | sin € 1/cos € I (Fessenkoff) | J (Schoenberg) 
° 0.00 I.00 I.00 I.00 
IO. 0.17 I.02 I.O1 1.04 
20. 0.34 1.06 1.03 ks 
30. 0.50 ts 1.07 I.34 
40 0.04 | £..2% I.13 1.05 
50. 0.77 1.56 r33 2.09 
60.. 0.87 2.00 1.36 2:32 
7°. 0.94 2.92 1.55 3.64 
8o. 0.98 5.76 1.85 5.00 
go I.00 oo 2.42 6.98 














The results are shown in Figure 3, where the values of J have been 
plotted against sin € as abscissa. The ratio Edge/Center is 2.4 for 
Fessenkofi’s formula, instead of 3.3 as computed in section 4 with 
¢.(a) = 0.5. The same ratio for Schoenberg’s formula is 7.0, instead 
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of 8.8, as computed in section 4 with ¢,(a) = sin? a/2. The conclu- 
sions of section 4 are not materially altered. 





° 0.2 0.4 0.6 0.8 1.0 Distance 
Int. T T T T from 
center 
6 





wa 
{ 




















There is an apparent discrepancy between the shape of the com- 
puted curves in Figure 3 and that of the observed curve in Figure 1. 
This is not serious, however. It is, of course, obvious from an ex- 
amination of the photographs that the curve does not reach a sharp 
point in the center of the core. The bottom of the curve must be 
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flat. I have attempted to get a microphotometer tracing from photo- 
graphs of the dark nebulae in Taurus; and while the result is not 
very satisfactory, the curve more nearly resembles Figure 3. I be 
lieve that more accurate observations will tend to confirm this re- 
sult. The steep edges shown in Figure 3 are probably cut down by 
the decreasing thickness of the nebula and by the failure to obey our 
assumption that R& is very large. 

6. While the majority of cores in dark nebulae follow the general 
outlines of the nebula as depicted from their obscuring effect upon 
background stars, not all dark nebulae have cores. This is very clear- 
ly demonstrated by the lanes designated as Barnard 8 in Plate VI. 
If our explanation of the phenomenon of dark cores is correct, we 
must conclude that a lane which fails to show a dark core is, as a 
whole, inclined to the plane at right angles to the line of sight. The 
fact that there is a conspicuous core at B11, on the left side of B8, 
suggests that B8 is an elongated cloud, the nearer edge of which 
merges with B11 while the visible surface of the cloud extends toward 
the west and away from us. 

7. In sections 4 and 5 we found that the observations favor a 
fairly large value of the spherical albedo, y, and a phase function 
which throws most of the scattered light in the direction of the 
incident beam. A normal phase function of the Lambert, Euler, 
or Lommel-Seeliger type can only be allowed if y is unreasonably 
large. It is difficult to estimate the uncertainty of this result. The 
amount and distribution of starlight may not be the same as that 
which is observed on the earth, although it is not probable that it 
will be materially different; the nebulae under consideration are 
not illuminated by nearby stars. The amount of the “terrestrial” 
background of the sky is rather uncertain, and so are the observa- 
tions with the Fabry photometer. Nevertheless, there seems to be 
a rather definite tendency in nebular investigations to give a value 
of y which is larger than that expected for ordinary spherical par- 
ticles of rock. The albedos of the asteroids average about o.1, 
while nebular particles give values ranging from about 0.5 to 1.0. 
This discrepancy does not necessarily mean that the material of 
the nebular particles is whiter than that of the asteroids. While 
in the latter we observe only the effects of diffuse scattering by 
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their surfaces, the particles of a nebula may be small enough to pro- 
duce appreciable diffraction effects. In accordance with Mecke," we 
may suppose that scattering by diffraction removes from the 
incident beam an amount equal to the cross-section of the particle. 
The total absorption coefficient is therefore 2Nzp?. Of the total 
absorbed light, each cubic centimeter of the medium scatters one- 
half as diffracted light. The other half is partly transformed into 
heat, and only the fraction yN 7p’ is scattered in accordance with 
the law of Lambert. If we consider only a single scattering process 
characterized by the phase function g(a), the spherical albedo y is 
an average between the ordinary value and unity—the latter repre- 
senting the 100 per cent efficiency of the diffraction process. 

8. One of the most interesting features of diffuse galactic nebulae, 
and one which suggests a definite phase effect, is the frequent oc- 
currence of bright rims which partly surround dark masses of nebu- 
lar material apparently projected against the structure of the bright 
nebula. The best known rim of this type is seen surrounding the 
Horse-head nebula, B33, from the northwest through the west 
and as far as the southwest. This rim does not belong to the under- 
lying structure of the bright nebula IC 434. 

Table 3 gives a list of the more conspicuous rims observed. It is 
very striking that, while many dark markings fail to show the rim 
effect, those that do show it always exhibit an exceedingly narrow 
crescent, often not more than o.o1 to 0.02 of the radius of the dark 
cloud in width. The transition along any such radius from dark to 
bright is exceedingly sudden. At the outer edge of the rim the bright- 
ness usually fades out more gradually. Wherever there is a definite 
symmetry in the underlying nebula, as in NGC 6523 or NGC 6514 
(Trifid nebula), the brightest portion of each rim is in the direction 
of the exciting star. This alone proves that the dark clouds with their 
rims form a part of the bright diffuse nebula upon which they are 
seen projected. Most of these dark objects seem to be somewhat 
closer to us than the exciting stars; but the difference in distance 
cannot be very great, because in most cases the dark clouds are not 
completely black but exhibit varying degrees of illumination. 

It is most surprising that the spectra of the nebulae exhibiting 


™ Ann. d. Phys., 65, 265, 1921. 
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bright rims are always of the emission type. Table 3, in which the 
nebular types are taken from Hubble, contains only one exception: 
IC 12744, ¢, is listed by Hubble as having a continuous spectrum, 


TABLE 3 
LIST OF BRIGHT RIMS* 





| | 
eee Fe) a 
ee | Exciting Sp. of | Sp. of | eee ae c ie 
Nebula i Chane | Star | Neb. | Description of Rims ource 
IC 434 ¢Orionis | Bo | E | Excellent rims on the northwest, | J. C. Duncan 
| o Orionis | Bo | | west, and southwest edges of 
| B33 (Horse-head nebula), of | Ap. J., 53, 392, 
the bay nebula near BD — | 1921 


2°1343 and of the small bay 
| nebula 6’ farther north. The 
entire structure of IC 434 may 

be considered as a rim of the 

large dark nebula to the east 


NGC 6523=M8 | Cluster Oes | E | Several fine rims pointing toward | Ross, Ap. J., 77, 
| Bo | exciting stars 246, 1933 
| | | Duncan, Ap. J., 
| | 51, 4, 1930 
| | | De Kerolyr, 
| | | L’ Astronomie 
| 48, 536, 1934 


NGC 6611 =M16. .| Cluster | Oes 1 E | Several fairly bright rims point- | Duncan, Ap. J., 
| | | ing toward brightest stars in | 51, 7, 1920 
| | cluster 
| 
NGC 6514=Mz20. .| —23°13804 | Oes | E | Several excellent rims in the Duncan, ibid 
| | southern part of the Trifid | 57, 142, 1923 
| | nebula pointing toward the | 
| exciting star | 
| | 
NGC 10977 ...| ef Orionis Br | C+E | Fairly well-defined bright border | Duncan, ibid 
| around dark bay north of | _ p. 138 
| bright nebula; points toward | 
| c’ Orionis 
NGC 7000 | a Cygni? A3p | C+E | Several bright rims pointing | Duncan,  ibid., 
| | } west and north p. 147 
Nebula preceding| | | 
NGC 7000... | a Cygni? A3p | E | Several faint rims surrounding | Duncan,  ibid., 
| dark nebulae | 63, 122, 1926 
| | 
IC 1274, b — 23°13907 Bz i North rim of Bor | Barnard, Aflas, 
| —23°139908 Br c | No. 30 
| | | | 
IC 1306... | +56°2617 Oes je Excellent border surrounding | Barnard,  ibid., 
| dark nebula covering stars | No. 49 
| +56°2598, +56°2604. Points 
| toward +56°2617, if rim is | 
| | considered to belong to dark | 
| (and not to bright) nebu- | 
| Josity 
| | 
NGC 1976 a Trapezium Oes5 | E | Several confused rims Ritchey, Yerkes 
Bo | Pub., 1, 380, 
| | 1903 
| 








* The original plates show numerous fine rims which are lost in the published reproductions. 


but this was merely inferred by him from the spectrum, Br, of the 
exciting stars. The nebular spectrum proper has not been observed. 
It is quite possible that the spectral types of the stars are slightly 
in error and that the nebular spectrum is of the emission type. All 
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other nebulae listed in the table, except IC 1396, have been observed 
spectroscopically to contain emission lines. IC 1396 is excited by the 
Oes star +56°2617, and there can be no doubt that its spectrum is 
also of the emission type. 

We conclude that the rims are strongest when e = go’, and are 
absent when 7 = go’. It is also obvious that we are not dealing with 
a normal phase effect, since the probability of observing so many 
objects at i Y 180° is too remote. Rather should we conclude from 
the symmetry of the structure of such nebulae as NGC 6523 that in 
most. cases 7 © go’. It is perhaps not an accident that rims are most 
often observed at rather large angular distances from the exciting 
stars. 

If we retain the theory of scattering in spite of the fact that the 
spectra of these nebulae show emission lines, we are confronted with 
the fact that the illumination of an opaque spherical cloud of par- 
ticles, according to the Lommel-Seeliger law, does not lead to very 
narrow rims except for 7 near 180°.’ The same is, of course, true in 
the case of Lambert’s law. 

Perhaps we can find an explanation of the rims if we return to the 
derivation of the original Lommel-Seeliger formula. Following the 
original procedure of Seeliger, we shall assume that the absorption 
coefficient for the incident beam is k,, while that for the scattered 
beam is k,. Then: 

cos 7 


T=uL ——_ o>—R(k: sec it+ka sec €) ; 
uLe(a) k, cose + k; cos 1 [x (3) 





Let us now assume that k, is very much smaller than k,. Then, if 
R(k, sec i + k, sec €) is very large, we may write 


cos 1 





t=? Le(a) 
ky k, . 
cos e+. cos2 
k, 
If k,/k, = 0, we have I « (1/cos €), which is exactly the effect ob- 
served in the bright rims. If &,/k, is very small, but not equal to 


zero, the rim will be less pronounced: for € = go’, 
ky 
I=C-. 
k, 
12 Tbid., p. 74. 
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If k,/k, = 100, the rim is already very greatly enhanced over the 
one which is shown in Figure 3 for Schoenberg’s formula. The third 
column of Table 2 gives the reciprocals of cos e€ which would apply 
to the case k, = o. 

Our explanation of the bright rims finds its most interesting con- 
firmation in the consequences of Zanstra’s theory of nebular radia- 
tion. It is now generally recognized that the absorption coefficient 
of a nebula for ultraviolet starlight beyond the Lyman limit is very 
large. After this radiation has been converted into visible radiation 
consisting largely of the Balmer lines and of the stronger forbidden 
nebular radiations, the absorption coefficient is greatly reduced® 
and the observed light passes relatively freely out of the nebula. 
Dr. Chandrasekhar has called my attention to the fact that it is not 
legitimate to put k, = o. This is, of course, also shown by the ob- 
servations: thus, B33 is not transparent to the visible light of IC 434 
beneath it. But it is probably correct to assume k, < k,. It is en- 
tirely possible that the ratio k./k, is not the same for all emission 
lines. It may be smaller for the forbidden nebular lines than for the 
Balmer lines. This suggests some interesting observational tests, 
which will, however, be difficult because of the faintness of the neb- 
ulae. 

The physical interpretation of the phenomenon is rather simple. 
Since k, is large, the incident starlight penetrates to a small depth 
into the nebula. But k, is small; hence we can see into much deeper 
layers of the nebula, provided the layers are excited to luminosity 
by the incident radiation. For e€ = o our ability to penetrate to 
great depth yields us nothing because the exciting light, even for 
small values of i, fails to reach the deeper layers. For € = go’, how- 
ever, we look lengthwise through a narrow edge, all of which is 
reached by the exciting radiation and through which we can see to 
great depths because &, is small. 

It is interesting to realize that from the photometric theory alone 
we could have concluded that k, « k, and that the visible radiation 
of the nebula is of an entirely different nature than the incident 
light from the star. 

We still have to explain the gradual diminution in surface bright- 


3 Rosseland, Theoretical Astrophysics, p. 324, 1936. 
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ness of the rims as we pass along a radius away from the center of the 
dark cloud. This effect is probably due to the fact that at the ex- 
treme edge of a cloud, R(k, sec i + k, sec €) in (3) can no longer 
be considered large. If it is small, we may expand the exponential" 
and use only the first term: 

IT = uRL¢(a) a 

COS € 

But this expression holds only for a nearly transparent plane-parallel 
nebula, in which R/cos ¢ is the distance along the line of sight from 
the nearer to the farther surface. The relation between J and 1/cos € 
is due to the increased path within the nebula, if observed slantwise. 
In a spherical nebula the factor R/cos € should be replaced by 
2R, cos €, where R, is the radius. Accordingly, J decreases toward 
the edge of a transparent nebula with cos e. This may even be accen- 
tuated by a decrease in nebular density toward the edge—a phe- 
nomenon accounted for in our formula by the factor » which con- 
tains V, the number of particles per unit volume. 

g. Two very remarkable bright nebulae appear in one of the most 
opaque parts of the dark cloud in Taurus. They are shown as Bro 
(a = ghr1™o, 6 = +2820) and Bri4 (a = 4532™3, 6 = +2595 
in Barnard’s Allas of Selected Regions of the Milky Way. Both nebu- 
lae are also shown in Plates IV-V. Bro is illuminated by a star of 
photographic magnitude about 15. The radius of the nebula is 
approximately 5’. It is evident that these values are not in agree- 
ment with Hubble’s expression’s 


m, + 5loga=11. 


In order to satisfy the formula the star should have a magnitude of 
7.5. Since the nebula exhibits little structure, it is not probable that 
the star is reduced in brightness by intervening clouds. The dim- 
ming of the starlight probably takes place in the nebula itself. 
B14 has no star brighter than photographic magnitude 17. The 
radius a = 1/3. Hence to satisfy Hubble’s formula we should have 


™ Schoenberg, op. cit., p. 36. 


5 Struve and Story, Ap.J., 84, 208, 1936. 
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m, = 10.5. Evidently in this case the star is also dimmed in the 
nebula by 6 or 7 mag. 

In both of these objects we are probably dealing with exception- 
ally dense and thick nebulae. The stars are so deeply imbedded in 
nebular material that only a small fraction of their light manages 
to penetrate. If 7 is the optical thickness, the fraction e~’ of the 
original starlight is observed. The diffuse light of the nebula, on the 
other hand, has a better chance to survive. If the nebular particles 
were perfectly white, so that y = 1, the total light of the nebula 
would obviously be 1 — e~’. If y is not equal to 1, the problem is 
more complicated. The analogous case of a planetary nebula has 
been solved by Chandrasekhar."© The problem of dense reflection 
nebulae is now under investigation by Dr. Henyey. 

YERKES OBSERVATORY 
February 19, 1937 


6 Zs. f. Ap., 9, 266, 1935. 


A PHOTOELECTRIC STUDY OF THE LIGHT 
FROM THE NIGHT SKY* 


C. T. ELVEY AND F. E. ROACH 


ABSTRACT 


Observations of the light from the night sky have been made with a recording 
photoelectric photometer. Records covering the entire hemisphere were made in about 
one hour. Observations were obtained at different hours of the night and at different 
seasons of the year, in order to determine the distribution of the various sources of 
radiation coming from the night sky. 

Observations of certain regions of the sky were reduced by analytical methods to 
evaluate the unknowns. All of the tracings, however, were reduced by a graphical 
method, in which all radiations having similar distributions over the sky were grouped 
with respect to their co-ordinate systems. 

Isophotes have been constructed for the zodiacal light, arranging the observations 
into four groups. There is a seasonal variation in the zodiacal light which might be in- 
terpreted as resulting from a cloud of meteoric particles which are in conjunction with 
the sun during the late winter months. 

An analysis of the light associated with the galactic system shows that after the 
effect of the stars has been removed there is an excess of light, which we shall call the 
galactic light, and which is probably produced by the scattering of starlight by in- 
terstellar matter. It also includes the galactic nebulae. The mean amount of galactic 
light at latitude o° is equal to 57 stars of the tenth photographic magnitude, or is 
equivalent to 5.6 mag. per square degree. 

Applying a theoretical discussion of scattering by Struve, it is found that the spheri- 
cal albedo is probably large and that in the galactic plane the medium is opaque 


The light from the night sky is composed of radiations coming 
from many sources, which we may enumerate as follows: direct 
starlight, scattered starlight, light originating in the upper atmos- 
phere, which is termed the “permanent aurora,”’ the polar aurora, 
the zodiacal light, light from galactic sources, which is probably 
starlight scattered by dust particles in interstellar space, and light 
scattered by our atmosphere from all these sources, as well as light 
scattered from terrestrial sources. 

The distribution of the component radiations over the sky varies 
with each source, and consequently the photometric study of the 
light from the night sky offers a rather complex problem. The in- 
tensity is a function of so many variables that an analytical solution 
for the various components is impossible. However, with the aid of 
approximations and with observations obtained under certain con- 

* Contributions from the McDonald Observatory, University of Texas, No. 3. 
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ditions it becomes possible to separate the various distributions 
and to determine them with some certainty. 

The problem of the light from the night sky may be attacked in 
three ways, namely, by spectroscopy, by studies of polarization, 
and by photometry. 

The spectroscopic results have been summarized recently by 
Georges Déjardin.t The emission spectrum of the night sky has 
received the most attention, since the feeble radiations are confined 
to narrow bands of wave-lengths and consequently are photographed 
much easier than is a continuous spectrum of the same total inten- 
sity. Déjardin discusses three regions of the spectrum—the visual, 
the photographic, and the ultraviolet. In the visual region the most 
prominent feature of the spectrum of the night sky is the green 
auroral line. In the photographic region, the one in which our work 
has been done (since our photometer contains a potassium photo- 
electric cell), a large number of radiations have been identified, 
mostly from the molecule of nitrogen. The most intense radiations 
are the Vegard-Kaplan bands of nitrogen, which have an excitation 
of 6.1 volts. Next in intensity are the bands from the first and second 
positive systems of nitrogen, with excitations of 7.4 and 11.0 volts, 
respectively. The fainter radiations are identified with those of the 
negative bands of nitrogen, the excitation potential being 19.6 
volts. The relative intensities of these radiations from the night 
sky are exactly the reverse of those observed in the polar aurora. 
In the ultraviolet region of the spectrum of the night sky, in addi- 
tion to the emissions, there is a continuous spectrum with absorption 
lines. 

The foregoing studies are from spectrograms taken in regions of 
the sky where the zodiacal light or the Milky Way are not bright. 
The spectrum of the zodiacal light has been observed by E. A. 
Fath,? who found it to be like that of the sun, showing that it is 
scattered sunlight. Fath has also observed the spectrum of the inte- 
grated light from the Milky Way and found it to be of the solar 
type. 

J. Dufay? has applied polarization methods to the study of the 


™ Rev. Mod. Phys., 8, 1, 1936. 
2 Bull. Lick Obs., 5, 141, 1900. 3 Bull. L’Obs. Lyon, 10, No. 9, 1928. 
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light from the night sky. The light originating from the upper 
atmosphere by emissions from atoms and molecules does not show 
polarization; but, on the other hand, the light scattered by the 
meteoric particles producing the zodiacal light should be polarized. 
Dufay has measured the amount of the polarization both in the in- 
tense parts of the zodiacal light and in the regions of the sky away 
from the ecliptic, and has shown that a part of the light over the 
entire sky can be attributed to an extension of the zodiacal light 
over those regions. 

The third method of studying the light from the night sky, the 
photometric, has been used by P. J. Van Rhijn,4 who made a long 
series of observations at Mount Wilson in 1913. The instrument 
employed is one described by L. Yntema.’ The observations were 
made visually by comparing the background of the sky with a screen 
illuminated by a known source of light. Calibrations were obtained 
from observations of stars, with the aid of an artificial star. The 
observations were made in circles parallel to the horizon, and each 
series was compared with the light from the sky at the North Pole. 
This comparison was necessary since there are rather large and ir- 
regular variations in the amount of light from the North Pole, 
presumably caused by the variation in the intensity of the perma- 
nent aurora. 

An inspection of Van Rhijn’s table of observations shows that 
about eight points of the sky can be observed per hour; thus, to 
cover the entire sky with sufficient detail requires many hours of 
work. With the changing intensity of the auroral light it is quite 
possible that a considerable amount of uncertainty will enter into 
the observational data. 

J. Dufay? has made extensive photometric investigations of the 
light from the night sky, particularly that coming from the region 
of the North Pole. He has used both a visual photometer, comparing 
the light from the sky with an illuminated screen, and a photographic 
method described by Charles Fabry® in 1910. The Fabry photome- 
ter is excellent for detailed studies of regions of the sky and gives 
precise results. However, in conducting an extensive survey of the 

4 Pub. Astr. Lab. Groningen, No. 31, 1921. 

5 Ibid., No. 22, 1909. ®°C.R., 150, 272, 1910. 
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sky, it has the same disadvantage as the visual observations, name- 
ly, it takes from 5 to 8 minutes to obtain a single observation. 

In the photometric analysis of the light from the night sky we may 
take the brightness of a given region of the sky to be the sum of the 
intensities of the individual radiations as mentioned above: 

S . dete. 2 x 


Bo = = + Ts + 


- Ss Sc : T Sc ‘ 
F(z) F(z) F(z) F(z) + SX c( ) + X Sc(z) (1) 


where all of the intensities are expressed as numbers of stars of the 
tenth photographic (or photoelectric) magnitude per square de- 
gree, and where B, is the observed brightness of a given region of 
the sky, S is the intensity of the direct starlight, A is the intensity 
of the permanent aurora at the zenith, ¢(z) is the law of variation 
of the permanent aurora with zenith distance, F(z) is the extinction 
factor, ZL is the zodiacal light at the observed region, X is the un- 
known galactic radiation, Sc(z) is the law of variation of scattering 
with zenith distance, and TJ is the total auroral, zodiacal, and 
galactic light over the hemisphere. The general solution of equa- 
tion (1) is readily seen to be impossible. However, by a reason- 
able choice of assumptions and observing conditions the unknowns 
may be determined. 

The contribution to the light from the night sky by the stars can 
be determined by star counts. However, such data are not available 
for all regions of the sky, and for the present we shall use the mean 
star counts as summarized by K. F. Bottlinger.’ The deviations 
from the mean number of stars for regions of the sky at galactic 
latitudes greater than, say, 30°—40° will be small, and for the pur- 
poses of the present investigation we may neglect them. 

The determination of the extinction factor, F(z), has been made 
with the photoelectric photometer from observations of bright 
stars. For zenith distances less than 75° the factor varies as sec 2, 
but at larger zenith distances the extinction deviates from the 
secant law. Our observations agree quite well with Bemporad’s 
extinction theory.*® 

The intensity of the auroral light at the zenith, corrected for ex- 
tinction, is no doubt a function of the physical state of the upper 

7Zs.f. Ap., 4, 370, 1932: 8 Handbuch d. Ap., 2, Part 1. 
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atmosphere. There is little doubt that it varies with time, as ob- 
servations of the integrated light from a region around the North 
Pole by Lord Rayleigh? and others show diurnal, seasonal, and long- 
period variations. Since it was necessary to make a series of observa- 
tions over several hours in order to determine the law of variation 
of auroral light with zenith distance, we have assumed that a cor- 
rection can be made for the diurnal variation from simultaneous 
observations of the light from a region of constant zenith distance, 
namely, the North Pole. 

Van Rhijn‘ has developed an expression for the variation of the 
auroral light with zenith distance, the function g(z) in equation (1). 
He assumes that the light originates in a layer rather high in the 
atmosphere, that the layer is uniform over the region from which 
the observations are made, and that the light received from this 
layer is proportional to the thickness of the layer in the line of vision. 
The height to the layer is called 4, which is expressed in units of the 
radius of the earth and is small in comparison with the latter. The 
thickness, €, is assumed to be small enough so that in the develop- 
ment of the function the squared terms in ¢ are negligible. The re- 
sulting expression for the function is 


gia) = : (2) 


sin? 2 


Vv? 7 (1 +h)? 


This function reduces to sec z when the value of / is taken as zero. 

It is seen that the intensity of the auroral light, Ag(s), increases 
toward the horizon, the rate of increase being greater for smaller 
values of h, the height to the auroral layer. Since the extinction is 
produced by the atmospheric layers between us and the auroral 
layer, its variation with zenith distance will be more rapid than the 
increase in the thickness of the auroral layer; consequently there will 
be a tendency for the extinction to overcome the increase in the in- 
tensity of the auroral light toward the horizon. This produces a 
maximum of intensity of the auroral light a few degrees above the 
horizon, which we have noticed visually many times. 

The zodiacal light is a function of the difference in longitude be- 


9 Proc. R. Soc., Ser. A, 1§1, 22, 1935. 
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tween the observed point and the sun, and of the latitude of the 
point. Numerous observations by visual observers, and our own 
records with the sky photometer, show that the zodiacal light is 
not symmetrical with respect to the ecliptic, the axis being as 
much as 5° north of the ecliptic. The work of both Van Rhijn 
and Dufay show that the zodiacal light is not merely confined to the 
region of the ecliptic but extends over the entire sky. Van Rhijn 
attributes a rather large portion of the light of the night sky to 
zodiacal light, while the studies of Dufay using polarized light 
indicate a much smaller percentage. 

The quantity listed as unknown galactic radiations includes the 
diffuse galactic nebulae and any scattered light coming from inter- 
stellar matter. It is probable that this light is concentrated in low 
galactic latitudes and may be neglected at the higher ones. 

The scattered light from all sources may be classed together, since 
we use the same scattering function, namely, that determined by 
C. G. Abbott’? from observations made at Mount Wilson. How- 
ever, since we have a value for the total number of stars, we shall cor- 
rect our data for scattered starlight before proceeding with the 
analysis of the observations. In correcting for the scattered star- 
light, we make the assumption that the stars are distributed uni- 
formly over the hemisphere; and we have modified Abbott’s values 
of scattering in accordance with the wave-lengths of light used by 
us, namely, those of the photographic region. 

The total auroral, zodiacal, and unknown sources of light which 
are scattered are grouped together as one term, 7. Again, the 
assumption is made that they are distributed uniformly over the 
hemisphere. 

Under the assumptions we have made, there remain four un- 
knowns for which we must solve: A, the intensity of the auroral 
light at the zenith; ZL, the zodiacal light at the point in question; 
T, the total light (other than from stars); and h, the height to the 
auroral layer. Since the height to the auroral layer enters into the 
function g(z) in a rather complex manner, it is much easier to solve 
for it by approximations, thus leaving three unknowns. 


10 A.J., 28, 129, 1914. 











































PLATE VII 


THE PHOTOELECTRIC SKY PHOTOMETER OF THE MCDONALD OBSERVATORY 
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APPARATUS AND METHODS 

A photoelectric photometer of the recording type was constructed 
for the purpose of making a comparatively rapid survey of the sky, 
thus eliminating one of the objectionable features of the visual 
observations. A series of tracings covering the entire sky could be 
made in about one hour. The sky photometer is shown in Plate VII. 
It is mounted on a table rotating on the vertical support, thus mak- 
ing it possible to observe in all azimuths at a given altitude. The 
lens forming the image of the sky on the photo-electric cell has its 
axis in the horizontal plane. By means of a mirror placed at an 
angle of 45° to the axis, light may be reflected into the photometer 
from any altitude. The photometer tube, a large brass cylinder, 
contains the Kunz photoelectric cell at one end and the amplifying 
tube, Pliotron FP-54, of the General Electric Company, in a com- 
partment at the other end. The amplifier is of the balanced circuit 
type designed by Du Bridge and Brown."' The amplifying tube, 
with its grid resistor (4 X 10" ohms), is located in the photometer 
tube, as near to the photoelectric cell as is possible. The remaining 
parts of the circuit are in the metal box shown on the table and 
are connected with the photometer by a shielded cable. The high- 
potential batteries for the photoelectric cell and the storage bat- 
teries for the amplifier circuit are in the box immediately behind 
the control box. 

The output of the amplifier is measured with a Leeds and Northup 
high-sensitive galvanometer, the period being 1 second and the 
current sensitivity 5 X 10 ? amperes per millimeter at 1 meter. The 
deflections of the galvanometer are recorded with a beam of light 
from a motion-picture exciter lamp for sound reproduction which is 
reflected from the galvanometer mirror to the recorder mounted on 
the vertical support of the photometer. The recorder consists of 
a shield with a horizontal slit, and directly behind the shield is a 
holder containing bromide paper. The holder is connected by a 
pulley system to the photometer so that it is lowered past the slit 
as the photometer is rotated. The motion of the bromide paper is 
such as to make half a millimeter on the tracing equal to 1° of azi- 


™ Rev. Sc. Instr., 45 532, 1033- 
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muth. Mounted at one end of the slit is a small light which is con- 
nected with contacts on the photometer so that the light is flashed 
as the instrument passes the cardinal points, thus giving the orienta- 
tion. A check on the orientation is obtained by making tracings 
across Polaris and other bright stars. 

A typical tracing is shown in Figure 1, a reproduction of a run 
made at a zenith distance of 70° on the morning of November 9, 1934. 
The tracing shows the zodiacal light in the east, the Milky Way, and 
several stars of medium brightness. The zero tracing shown at the 
bottom of the figure is placed on each sheet, inasmuch as it shows a 
departure from a straight line. 
The procedure is to make a short 
zero tracing with the photometer 


. to the north, then open the shut- 

a ter and make the circuit of the 

~~ / sky, and close the shutter and 
reverse the run to make the com- 


plete zero tracing. The apparent- 
rel ww} ly Sinusoidal variation in the zero 

Viti saint tracing was found to be due to 
the insufficient magnetic shield- 
ing of the Pliotron tube from the 
earth’s magnetic field. Through- 
out any one night the zero tracings were reasonably constant; how- 
ever, from night to night there were slight differences in amplitude, 
although the shape of the curve was the same. 

The observations consist partly of tracings covering the entire 
hemisphere and partly of a series of observations of individual 
regions of the sky for the purpose of determining the variations of 
the auroral and scattered light with zenith distance. Observations 
were taken for the purpose of calibrating the instrument, as well as 
for the determination of the extinction factor. All observations 
were made long enough before or after the rising or setting of the 
sun or moon so that no scattered light from these sources would 
interfere with the results. The observing station was near Mount 
Locke in the Davis Mountains, at an altitude of 6200 feet; and the 
distance to the nearest source of terrestrial illumination was 10 


N E S W N 














Fic. 1.—A tracing made with photo- 
electric photometer on November 9g, 1934. 
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miles to Fort Davis, Texas, a small town without street lights and 
hidden at a lower altitude, behind the mountains. The nearest 
towns with street lights were: Alpine, 32 miles, Marfa, 26 miles, 
and Pecos, about 60 miles, all hidden from view by the mountains. 

The observations were made on the best photometric nights; 
and if any change was seen in the appearance of the sky, it was noted 
in the records. 

The normal procedure in securing a set of tracings of the sky was 
to observe a number of standard stars and regions of the sky before 
and after taking the set. The tracings were usually taken at zenith 
distances of 75°, 70°, 60°, 50°, 40°, and 30°. In general, it was our 
intention to obtain a series of observations for the evening and for 
the morning portions of the sky, each dark of the moon; but this 
was not always possible. Other series were taken at definite orienta- 
tions of the sky in order to obtain the zodiacal light or the Milky 
Way under better conditions. At sidereal time 12"40™ the Milky 
Way is practically on the horizon in our latitude, with the result 
that most of the entire set of tracings can be used to study the 
zodiacal light. For the study of the galactic light, tracings were 
taken when the Milky Way passed through the zenith near mid- 
night, when the zodiacal light is weak. At the same time the plane 
of the Milky Way is nearly vertical, and the tracings are thus cross- 
ing the Milky Way at right angles. The auroral light for each tracing 
is constant, and the variation of intensity is practically all due to 
variation in galactic light. 

To determine the extinction factor, observations were made on 
bright stars, a Lyrae, a Aquilae, a Virginis, and others, over a large 
range in zenith distance. Since the photometer received light from 
a large region of the sky around each star, it was necessary to obtain 
observations of the sky background on either side of the star. 
These background observations also gave data for studying the 
variation of the auroral light with zenith distance. 

Additional observations of given regions of the sky over a large 
range in zenith distance were made for use in solving equation (1). 
The regions chosen were those for which star counts could be had, 
such as certain selected areas at high galactic latitudes. Two of the 
principal regions used were areas 52 and 58. Also, observations were 
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made of the north pole of the ecliptic. Throughout a series of ob- 
servations, we observed the light from the North Pole at regular 
intervals in order to obtain a correction for the diurnal variations of 
the auroral light. 

An examination of equation (1) shows that the observations taken 
near the horizon will have by far the greatest weight, and will be the 
ones determining the accuracy of the solution. Also, the observa- 
tions made near the horizon will be subject to larger observational 
errors caused by small amounts of haze. For these reasons it seems 
best to combine observations of various regions of the sky into one 
least squares solution. This introduces the question of the variation 
in the amount of zodiacal light from place to place, as well as of the 
galactic light, when the observations are near the Milky Way. We 
have made corrections for these variations by reducing all observa- 
tions to one region of the sky by the following empirical method. 
We made a diagram of the sky brightness plotted against zenith 
distance, drawing a smooth curve through the observations for 
each region of the sky. We then interpolated for the brightness of 
that region at any given zenith distance—say that of the North 
Pole. Under our assumption, the auroral light will be the same for a 
given zenith distance, as will also the scattered light. The variations 
in these brightnesses from different regions of the sky can be inter- 
preted as being due to differences in zodiacal light and galactic light. 
From these data we can reduce all of the observations to one region 
of the sky, thus giving a larger number of points to be used in the 
least squares solution of equation (1). 

As mentioned above, the function ¢g(z) in equation (2) contains an 
unknown, h, the height to the auroral layer; but since it enters the 
function in a rather complex manner, it is necessary to solve for it 
by some other method than least squares. We have assumed values 
of h varying from 0.00 to 0.03. With each, a least squares solution 
was made for the other unknowns, from which we determined the 
residuals. The assumed values of / were then plotted against the 
sum of the squares of the residuals, and the point of minimum gives 
us the best value of h. 

We have made analytical solutions of the foregoing equation for 
data from two nights in which there was considerable difference in 
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extinction, principally to guide us in the graphical evaluation of 
the major portion of our data to be described below. The values of 
determined for these two nights are shown in Figure 2. The data 
for February 21, 1936, give a sharp minimum at a value of / equal 
to 0.008, or about 50 km; while the data from the night of October 2, 
1935, give a broad minimum at o.o1, or 65 km. 

These heights are some- 
what lower than the mini- 
mum height usually ob- 
served for the polar aurora, 
that is, 80 km. No doubt °%°°° 
we should expect the layer 
of the permanent aurora to 
be located in a different <5 00 
region of the atmosphere 
than the polar aurora, since 
there is a marked difference 
in the relative intensities of 
the radiations in the two 
spectra. It is also of inter- 
est to note that the ob- 30000 
served heights of the auro- Paes — agen pire 
ral layer are somewhat low- Fic. 2.—Graphical solutions for /, the height 
er than the height found to the auroral layer. Abscissae are the values of 
a en frequency af — niga ae a of evade ie “ the 

residuals. The dots are for observations made on 
meteors, which is also October 2, 1935; and the circles are for February 
about 80 km. ay 7907 

Since the radiation known as the permanent aurora must originate 
in a region where the atoms or molecules are in an excited state, it 
is tempting to assign them to the region of the Kennelly-Heavyside 
layer. However, measurements made of the virtual heights to this 
layer (or layers, as we know now) show them all to be greater than 
that we have found for the auroral layer. The ionosphere is com- 
posed of several layers, and the height of each may be determined 
by the reflection of radio signals. Measurements show that the 
lower, or E layer, is too—120 km in height, the F, layer, 200-230 km; 
and the F, region, about 240-500 km.” These meager data indicate 
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that the auroral layer is considerably beneath the layers which 
reflect radio waves. It is impossible for us to draw definite con- 
clusions from only two observations compared with the mean. It 
might be desirable to obtain simultaneous observations of the 
heights to the auroral layer and to the lower ionospheric layer. 

The details of the solutions for the two nights of observations will 
not be given. The evaluation of the unknowns, A, ZL, and 7, in 
the solution giving the best values of #4 is shown in Table 1. The 
agreement is closer than might be expected, except for 7, the total 
light. The large value of T for October 2 is no doubt due to our 
failure to allow for an increase in the numerical value of the scatter- 
ing function for a different value of the extinction. We have asumed 
that the variation of scattering with zenith distance is the same for 
each night; and consequently, the quantity TJ contains a factor 
which represents the variation in the scattering coefficient. 
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TABLE 1 
| | | 
Date | A | ZL (Pole) | T 
BEE : —| kid gs : 
Oct. 2. :2035).;,. | 133 57 | 1D. K IC 
Feb. 21, 1936. 136 | 45 








The analytical solution has been applied only to the observations 
of individual regions of the sky. The majority of our observations 
are tracings which recorded the brightness of the sky in a circle 
parallel to the horizon. The reduction of these observations by 
an analytical method would be rather tedious; and since we are 
especially interested in studying the light which may originate 
from scattering in the galaxy, we have used a comparatively rapid 
graphical method. We shall group as units the sources of light whose 
distributions over the sky are functions of the same co-ordinate 
systems, that is, the alt-azimuth, the ecliptic, and the galactic co- 
ordinates. 

We may best illustrate the method by following the reductions of 
an example. For this we shall take the observations of the evening 
of January 27, 1935, which afford a good example of the method we 
use in correcting for the aurora and for scattered light. Unfortunate- 
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ly, the tracings for this night were somewhat fogged by light during 
development and are not suitable for reproduction. The tracing in 
Figure 1, however, is typical; and we may refer to it. The abscissae 
of the tracings are azimuths, and the ordinates are deflections of 
the galvanometer. The lower tracing is that of the zero (no light 
on the cell), which, as mentioned before, is not a straight line, be- 
cause of the influence of the earth’s magnetic field. The rather 
sharp deflections in the zero tracing which are also present in the 
tracing for the sky are probably produced by ionization by cosmic 
rays in the air in the vicinity of the photoelectric cell and the ampli- 
fying tube. A. E. Whitford's found that it was necessary to evacu- 
ate the cell and the tube chamber in order to eliminate this un- 
steadiness. We decided against the evacuation of the apparatus 
because the equipment was in the field and because we were able to 
recognize the sporadic deflections on the tracings and to draw smooth 
curves to eliminate them. This is the first step in the reduction of a 
tracing. The next step is the elimination of the effects of the bright 
stars. A small celestial globe was set in the corresponding position 
of the sky at the time the observation was made, and with the aid 
of a spherical protractor and a knowledge of the size of the field 
of view of the photometer, about 9 square degrees, we could easily 
identify the deflections due to the brighter stars. All individual 
stars brighter than the third magnitude are eliminated with cer- 
tainty. 

After having drawn a smooth curve through the tracing, which 
represents the sky background, we next measured the tracing with 
the aid of a millimeter réseau, recording the deflection at every 5° 
of azimuth. The azimuth scale of the tracing was constructed so 
that each millimeter equals 2°. 

Having measured the tracing, we must convert the deflections of 
the galvanometer into known units of intensity, which are numbers 
of stars of the tenth magnitude per square degree. We do this with 
the aid of the bright stars used for calibration and extinction de- 
terminations. We did not have accurate photoelectric magnitudes 
for these stars, and the photographic magnitudes were not satis- 
factory, many of them having been determined by applying the 


3 Ap. J., 76, 213, 1932. 













































226 C. T. ELVEY AND F. E. ROACH 


mean color index of each spectral class to the visual magnitude. We 
have thought it best to apply the photoelectric colors of the indi- 
vidual stars to the visual magnitudes, which has given us fairly 
consistent results. For the photoelectric colors we have used the 
measures by K. F. Bottlinger™4 and W. Becker,’s and the visual mag- 
nitudes are those of the Harvard Revised Photometry. We have 
checked the resulting magnitudes with our observations by com- 
paring the standard stars with Polaris and applying a correction 
for the variation in light of the latter. The photoelectric magni- 
tudes are not accurate enough to be listed as such but are sufficient 
for the purpose of this investigation. 

The co-ordinates of the points which have been measured are all 
in the altitude-azimuth system. It is necessary to know the positions 
expressed in the other systems, and for this we have used graphical 
or mechanical methods of conversion. In changing from right as- 
cension and declination to galactic co-ordinates, we have used the 
tables computed by P. Emanuelli.”° 

The observed intensities have all been corrected for the intensity 
of the direct starlight by making use of the tables of mean numbers 
of stars for various galactic latitudes by K. F. Bottlinger.'? He has 
expressed his data in units which are numbers of tenth-magnitude 
stars per square degree. To make this correction it was necessary 
to take into account the extinction, which has been determined for 
each night from the observations of the standard stars. As men- 
tioned before, this correction for direct starlight has been applied 
for all galactic latitudes; but the mean data are rigorous for lati- 
tudes greater than 30° or 40°. 

A correction is applied, also, for the scattered starlight. The law 
of scattering with zenith distance is that determined by G. C. 
Abbott,” and the numerical value has been adjusted to the region 
of the spectrum to which the photoelectric cell is sensitive, namely, 
about 4500 A, assuming that the scattering is primarily according to 
Rayleigh’s law. This is probably very close to the truth, since we 
have observed only on the best photometric nights; and at this 


14 Veréf. Berlin-Babelsberg, 3, Part 4, 1923. 
5 [bid., 10, Part 3, 1933. 
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height above sea-level the sky is usually blue during the day, even 
near the horizon. As mentioned above, the additional assumption 
is made that the stars are distributed uniformly over the hemi- 
sphere; and we have taken for the total number of stars the value 
given by Seares, Van Rhijn, Joyner, and Richmond."* The correc- 
tion for scattered starlight rarely, if ever, amounted to more than 
4 or 5 per cent of the observed sky brightness. 

After corrections have been applied for the direct and the scat- 
tered starlight, the data from each tracing parallel to the horizon 
are composed of a constant amount of auroral light (the so-called 
“permanent aurora’’), a constant amount of scattered light, pos- 
sibly a variable amount of polar aurora, and variable amounts of 
zodiacal light and radiations associated with the galaxy. Since we 
are primarily interested in determining the amount and distribution 
of the radiations associated with the galaxy, we can group the two 
sources of illumination which are functions of zenith distance and 
solve for the two together. This is readily done by a graphical 
method, since we can make the assumption that all the galactic 
radiations are at latitudes, say, less than 30° or 40°; and for higher 
latitudes the correction has been made for the stars. We are also 
reasonably certain that in the regions away from the sun, the varia- 
tion of the zodiacal light over a short distance in longitude or lati- 
tude is rather small, and as a first approximation it may be con- 
sidered constant. A part of the tracings for each zenith distance 
made on a given night cross the regions just mentioned, which lie 
in high galactic latitudes far from the sun. Within this region there 
should be a relatively constant amount of zodiacal light for which, 
as a first approximation, we can use the value determined for the 
North Pole or for the pole of the ecliptic. 

The solution is best carried out by means of the diagram shown 
in Figure 3, which is a plot of the observed points on an equal area 
projection of the sphere. The co-ordinates are longitude and lati- 
tude. Every other observation has been omitted for convenience in 
making the illustration. The galactic equator is drawn, as are also 
the lines representing the limits of + 30° galactic latitude. The ob- 
served intensity for each point is also given. By referring to Figure 3 
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it will be seen that the region north of the ecliptic at a longitude of 
180° from the sun is far from the galactic plane, and that the vertical 
circles are roughly parallel to the circles of latitude. This will give 








Fic. 3.—Plot of observations for January 27, 1935, to illustrate the method of de- 

















termining the auroral and scattered light and the isophotes of the zodiacal light. } 
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Fic. 4.—The variation of the aurora plus scattered light with zenith distance for 
the observations shown in Fig. 3. 
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a number of points at different zenith distances, while the zodiacal 
light remains roughly the same. Assuming values of the zodiacal 
light for the first approximation, we are able to plot the points cor- 
rected for the zodiacal light in Figure 4 with the corrected intensities | 
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as ordinates and the zenith distances as abscissae. A smooth curve 
drawn through the points represents the variation of the sum of the 
permanent aurora plus scattered light. By correcting each set of 
observations on a circle parallel to the horizon by the respective 
amounts of auroral and scattered light, we obtain the values of the 
zodiacal light for the regions of the sky outside the limits of + 30° 
galactic latitude. These isophotes are shown by the continuous 
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Fic. 5.—Plot of observations from one tracing made on January 27, 1935, after 
correction for auroral plus scattered light. The continuous curve shows the variation 
of zodiacal light obtained from Fig. 3. The difference between the points and the curve 
gives the galactic light. 


lines in Figure 3. We can readily interpolate across the Milky Way 
for the intensities of the zodiacal light in that region. 

Our data for any given tracing parallel to the horizon, corrected 
for scattered light plus auroral light, may then be plotted as in 
Figure 5, with the intensities as ordinates and the azimuths as 
abscissae. A continuous line is drawn through the points represent- 
ing the intensity of the zodiacal light, obtained from Figure 3; and we 
can readily see that in crossing the Milky Way (azimuths 135° and 
320°) there is an excess of light. This excess will be called the galac- 
tic light; it is caused by the deviations of the actual numbers of 
stars at each point from the mean for the corresponding galactic 
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latitude, by the light which originates in galactic nebulae, and by 
light scattered by interstellar matter. 
THE POLAR AURORA 


The plot of the observations made for the purpose of determining 
the galactic light sometimes (Fig. 5) shows an excess in the region 
of the north point. The polar aurora is very seldom seen at the lati- 
tude of Mount Locke. We have suspected occasionally that the 
northern horizon was slightly brighter than other regions, but at no 
time during our series of observations have we seen anything that 
might be called an auroral arch. However, on several of our tracings, 
as just mentioned, we have noticed an excess of light in the north 
which cannot be explained. In particular, the nights of January 25 
and 27 and February 26 and 27 have shown a rather marked 
maximum around the north point for the tracings made at zenith 
distance 60° or greater. As a rule, the excess is larger for the tracing 
of greatest zenith distance and decreases with decreasing zenith 
distance. We have measured the point of maximumintensity on 
each of the tracings, and the mean is at azimuth 188° +7° (m.e.). 
The azimuth of the magnetic north pole is about 192°. 

THE ZODIACAL LIGHT 

By the method described above we have obtained the isophotes 
of the zodiacal light for each night for which tracings have been 
made. On January 9, 1935, the set of observations was obtained 
when the Milky Way was on the horizon. This gave a good picture 
of the zodiacal light for the hemisphere. We made use of this in- 
formation to reconstruct the regions of equal zodiacal light on other 
nights, in order to obtain the relationship of the auroral and scat- 
tered lights with zenith distance. It was soon evident that the zodi- 
acal light was not constant with the seasons, and this night could 
be used as a guide for only a short interval. Other nights were chosen 
as we worked away from this date. It might be remarked that if our 
assumptions about the regions of equal zodiacal light were not ex- 
actly correct, we would obtain a scattering of the points. An em- 
pirical correction could be made if one tracing gave results con- 
sistently too high or too low with respect to the others. 

The observations have been arranged in four groups. Figure 6 
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shows the data for six nights in November and December, 1934. 
Figure 7 is for six nights in January, February, and the first week 
in March, 1935. Figure 8 is for five nights in April and May, 1935, 
while Figure 9 contains the results for three nights in June and 


July, 1935. 
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Fic. 7.—Isophotes for the zodiacal light for six nights in January, February, and 


March, 1935. 


An inspection of Figures 5-8 shows the seasonal changes men- 
tioned above. In November and December the evening zodiacal 
light was rather inconspicuous and the intensity increased toward 
the sun rather slowly, while the gradient of intensity for the morn- 
ing zodiacal light was steeper. The isophotes for January and 
February show gradients of intensity with longitude of about the 
same order for the evening as for the morning light, with the morn- 
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ing light having a slightly greater intensity. In the figure for April 
and May we see that the evening zodiacal light has increased greatly, 
while that for the morning has decreased. The gradient has become 
quite steep for the evening light and distinctly weak for the morn- 
ing light. The tracings for the months of June and July have the 
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Fic. 8.—Isophotes for the zodiacal light for five nights in April and May, 1936 
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Fic. 9.—Isophotes for the zodiacal light for three nights in June and July, 1936 


zodiacal light only for the morning. It is seen, however, that there is 
an increase in the gradient. 

Another way to illustrate the changes in the zodiacal light is to 
take sections across the zodiacal light at longitudes from the sun, 
say, equal to 60° (evening zodiacal light) and to 300° (morning 
zodiacal light). These sections are shown in Figure 10, which has as 
abscissae the latitudes and as ordinates the intensity of the zodiacal 
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light. For the evening zodiacal light the greatest intensity is for the 
April-May series, while the least intense is the January-February 
series. The morning zodiacal light shows the reverse of the fore- 
going. 

From the data we might conclude that there is a cloud of particles 
of greater than the average density in such a position that it is in con- 
junction with the sun during the late winter months. However, the 
data are too meager to justify an accurate determination. An accu- 
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Fic. 10.—Sections across the evening and morning zodiacal light, at \ — A, = 60° 


and 300°, respectively. The heavy continuous line represents the data for April and 
May. The heavy dashed line is for January and February. The narrow continuous 
line is for November and December; and the dotted line is for June and July. 


rate set of the isophotes of the zodiacal light for each dark of the 
moon during the winter and spring months would undoubtedly give 
interesting information concerning the size and the parallax of the 
cloud. 

THE GALACTIC LIGHT 


In the preceding sections we have shown how to remove those 
components of the light from the night sky which are functions of 
the azimuthal and the ecliptical co-ordinate systems. The residuals 
are functions of the galactic co-ordinates. Since each of our ob- 
served points represents the light coming from an area of 9 square 
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degrees, there has been sufficient integration to obliterate the finer 
details of the Milky Way, and we have decided, for convenience in 
the construction of the isophotes of the galactic light, to average all 
points falling within a square of 5° on a side. This has been done 
in Figure 11, in which the isophotes were drawn from the mean 
values for the excess light plotted in the centers of every 25 square 
degrees, using the galactic co-ordinates. The first thing noticed is 
the lack of deficiencies. If we were dealing only with the deviations 
from the mean numbers of stars, then the plot should show nearly as 
many deficient areas as areas of excess. Before discussing the ob- 





Fic. 11.—Isophotes of the excess light associated with the galactic system. This 
represents intensities due to deviations from the mean number of stars, to galactic 
nebulae, and to starlight scattered by interstellar matter. 


servations, we must consider two assumptions which have been 
made in the reductions of the observations, namely, (1) that the 
light being scattered by our atmosphere is uniformly distributed 
over the hemisphere, and (2) that mean star counts were applied 
at low galactic latitudes. 

Since a large number of observations have been obtained with the 
Milky Way in different positions in the sky, we have investigated the 
dependence of the excess light for a given region upon the zenith dis- 
tance of the object during the observation. As the scattering is pro- 
portional to the air mass, which in turn is proportional to sec z, the in- 
tensities have been plotted against sec z. This should give an approxi- 
mately straight line which can be extrapolated to sec z = 0, to give 
the correct value. The intensities of the various regions have been 





THE LIGHT FROM THE NIGHT SKY 


iS) 
os) 
wn 


adjusted so that the mean intensity is too. Figure 12 shows such a 
diagram in which have been combined all data covering the Cygnus 


region of the Milky Way. We 
see from it that we should re- 
duce our observed 
light by 35 per cent. This 
amount seems rather large. 


galactic 


However, if one considers that 
it includes scattered light from 
the stars in the Milky Way, 
as well as the galactic light, 
which is nearly as intense, the 
value seems more reasonable. 
The data shown in Figure 11 
have been corrected by the 
foregoing amount. 


_ —— r . ——y a o 





ral 
«| 140 


» — 4[20 


0.00 1.00 2.00 3.00 4.00 


Fic. 12.—Variation with sec z of the excess 
light in the galactic system. Extrapolation to 
sec s = o shows that a correction must be ap- 
plied to the observed excess light. This cor- 
rection has been applied to the data plotted 


in Fig. 11. 


Figure 11 shows the excess light to be strongest in the region of 
the Sagittarius cloud, longitude 300°, and to be faintest—treally a 
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in the region of Auriga, longitude 120°. 
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Since there are comparatively few individual star counts in the 
Milky Way, it is better to study the mean intensity of the excess 


light with respect to latitude. 


The mean distribution of the excess 


light, according to galactic latitude, is shown in Figure 13 and in 
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Table 2, Groups 1-8. The small circles and the dotted curve show 
the mean observed excess of galactic light over the mean number of 
stars, for the whole of our data. This, however, is not a fair test of 
the excess light, since our data cover the regions where the stars 
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Fic. 13.—Plot of mean galactic light against galactic latitude. The dashed curve 
represents all the data. The continuous curve is for data from the same regions for which 
the mean star counts have been made. The individual points are for areas for which 
actual star counts have been used, and represent true galactic light. 


are most numerous, namely, the Sagittarius cloud, while the data 
for the mean star counts do not include this region. A further com- 
parison was made by grouping our data in such a way as to include 
only those which fall into the regions or groups of selected areas used 
in the determination of the mean number of stars."* The mean 
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excess is then shown by the continuous line in Figure 13, and by 
Groups I-IV in Table 2. 

For certain selected regions we have individual star counts; con- 
sequently we can correct our excess for the observed deviation of 
the number of stars from the mean. Bottlinger’’ gives the total 
brightness of the direct starlight for the Sagittarius cloud, for the 
Cygnus cloud, both bright and dark, for a region in Ophiuchus, and 
for the Scutum cloud. Having corrected our data for these particular 
regions, for the deviation of the number of stars from the mean, we 
obtain the excess light in each region as shown in the diagram. 

We were able to obtain data on the total amount of starlight from 
fourteen of the selected areas, making use of the information given 
in the Bergedorfer S pektral-Durchmusterung. Corrections to the ex- 
cess in these regions for the deviations of the actual number of stars 
from the mean give the true excess. These have been plotted in 
Figure 12 by a dot with the number of the selected area beside it. 
All of the selected areas show an excess, with the exception of area 
10, which has a small deficiency. This area is in Camelopardalis, 
just off the edge of the Milky Way where it is of the lowest density. 
In the chart of the isophotes of the galactic light this point falls in 
the region of deficiencies at longitude +94° and +11°. 

The data are shown also in Table 3. All the regions show an 
excess after all known sources of light except the galactic nebulae 
have been removed. The rather outstanding result that the dark 
regions show a greater excess than the regions of larger star density 
may be partly illusory since all points falling within a square of 5° 
have been averaged. This will obliterate much of the detail in the 
structure of the Milky Way. However, this smoothing of the ob- 
servations is hardly enough to explain all the excess in some of the 
larger dark areas. Also, we should expect the dark areas to be some- 
what luminous since they are clouds of cosmic particles which are 
obstructing the light of the stars behind them and are illuminated 
by the stars, much in the same way that a thin cloud in the night 
sky is illuminated by sky radiations rendering it no darker than the 
sky, although it obstructs a rather large amount of starlight. 

Observations by Struve and Elvey’® of some of Barnard’s dark 


19 Tbid., 83, 162, 1930. 
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nebulae gave results which indicated considerable scattering of 
starlight by the material in interstellar space. Struve?? has dis- 
cussed the Rayleigh scattering of starlight by interstellar matter. 
Assuming that the number of particles per cubic centimeter is that 





TABLE 3* 








IGalactic Light 


Galactic Corrected for 





oe | Latitude _|the Actual Num 
| ber of Stars 
Region around Ophiuchus 
cloud. -4 | 20 
| | 
Ophiuchus cloud, dark + 2 77 
| 
Cygnus cloud, bright + 6 78 
Cygnus cloud, dark.... fe) 116 
Scutum cloud..... | —4 | 32 
Sagittarius cloud | — 7 90 
Selected Area | 
No. I +27 18 
No. 2 +12 40 
“9 os | jie 10 
ee 31 8 
No. 6... +35 I 
No.7... +20 28 
9 Ser : 3 36 
No. 9. 2 23 
No. 10 +11 — 14 
BOLE... +25 4 
NO. 12. < . | +40 14 
NOE? . .. 2 seca +18 35 
No. ee : | + 6 25 
ae ay — 2 60 











* The values in the last column are expressed in terms of the equiv 
alent numbers of stars of magnitude to per square degree. 
given by Gleissberg,7"” he obtained a rough value for the surface 
brightness of the galactic light which is somewhat greater than the 
observed value. 
Struve’s more recent results” can now be applied to the discussion 
of our observations, and we are indebted to him for the following 


20 Thid., 77, 153, 1933- 
1 4.N., 246, 329, 1932. 2 A». J., 85, 1937. 
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treatment. We may consider that each volume element in inter- 
stellar space receives, on the average, the same amount of starlight 
that an equal element in the solar system would receive from the 
entire celestial sphere. In other words, we shall assume that the 
surface brightness of the starlit sky as seen from each volume ele- 
ment is uniform and constant, equal to the surface brightness pro- 
duced by all stars, as seen from the earth. Let this amount, L, 
be given in terms of the number of tenth-magnitude stars per 


square degree. Thus 


L, = 1953 Ld = 577 Stars of magnitude 1.o. 
T 
Obviously, 
Lo j 
L= = 56 stars of magnitude 10.0. 
41,253 
If only first-order reflections are considered, and if y is the spherical 
albedo of each particle, p is its radius, V is the number of particles 
per cubic centimeter, and R is the extent of the diffuse medium in 
the line of sight, the surface brightness of the galactic light is given 
by 
I, = Lylt — e~Xtrk] , 


provided the phase function g(a) = 0.5. J; is expressed in num- 
bers of tenth-magnitude stars per square degree. 
The observations give J = 57. Hence 


y[1 aa e~NxoR| =I, 


This can only be satisfied if 
Nrp’Rwo and y¥rSI. 


In other words, the medium extends to very large distances R, so 
that the transmission, measured by e~*"”*, is effectively zero. 
Furthermore, the albedo is unexpectedly large. It should be empha- 
sized here that the observational value, J = 57, is not very precise, 
and that y measures not only the efficiency of ordinary scattering 
but also the fluorescence which give rise to emission nebulae. 
Nevertheless, it seems desirable to extend the formula for J so that 
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it will include the second-order reflections. For this purpose we shall 
assume that each volume element in space, dv, is surrounded by an 
infinitely large sphere consisting of volume elements dv’, each of 
which receives the same amount of light, (41,253/47){Ldw, from 
the starlit sky. If we again assume that the phase function g(a) = 
0.5, then the element dv’ contributes the following amount of light 
to dv: 


e—kA 
dq = 0.5m f Lda Ae dv dv’ , 


where A is the distance between dv and dv’. The other designations 
are the same as those used by Struve.” Of this amount, dv reflects 


toward the observer, at € = 0, 


e—kA—kr , 
dq: = 0.25u? f Ldu A? dvdv'. 


Since 
dv = deadr, 


dv’ = A’ sin @dAdgdé , 
we have 


QO = 0.25p? do f Ldw {Sf e—*4—krdy sin Od0dAd¢ . 


We integrate over the entire sphere, and from A = 0 to A= « 
and fromr = otor = «. Struve” has shown that 


p=N 7 Tp , 
21 
k = Nap 


Hence 


I, = 0.25 . (4nr)/*L = y*L. 


Adding J, to J,, we have 
l=1,+ 1, = Ly(1 +7). 


This would give y = 0.6, if we again take 7/L = 1. 
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It would not be profitable to continue this speculation by building 
up a power series in y, for we have already made the untenable 
assumption that the diffuse matter extends to infinity in all direc- 
tions from dv. Since in reality the scattering layer is thin, we cannot 
carry out the integration to A = oo in all directions. It is not pos- 
sible to determine y without more information concerning the ex- 
tent of the medium, and concerning the values of p and N, which 
become important when the upper limit of integration in A is finite. 
The only conclusion we can derive from our present data is that 7 
is not small enough to permit us to ignore completely the reflections 
of higher orders than the first, and that in the galactic plane Nap?R 
is large enough to consider the medium as opaque. This agrees well 
with our current conception of the absorbing layer which gives rise 
to the zone of avoidance of extragalactic nebulae.” 


We are indebted to our colleagues of the Yerkes Observatory for 
many helpful suggestions and for their assistance, and to Dr. F. H. 
Seares for information concerning star counts used in this paper. 
Mr. L. Henyey assisted in the construction of the photometer, and 
Mr. Theodor Immega has given valuable service in obtaining the 
observations. We are indebted to the Allan-Bradley Company, of 
Milwaukee, Wisconsin, for one of the high resistances used in the 
photometer. 

McDONALD OBSERVATORY 
UNIVERSITY OF TEXAS 
February 25, 1937 

3 After this paper was written we received fasc. 19 of Vol. 1 of the Publications de 
V’Observatoire de Lyon, in which Wang Shih-Ky derives a theoretical value for the 
diffuse galactic light of about one-half of L. He has assumed Rayleigh scattering and 
has obtained the scattering coefficient from the observed value of the interstellar 
absorption. Considering the fact that his upper limit for J is equal to L, we can con- 
clude that the theoretical data agree, as to order of magnitude, with the observations. 


Note.—In equation 1, page 216, the fifth member on the right-hand side of 
the equation should be Sr X Sc(z) instead of S X Sc(z) since in this case the 
quantity is the total number of stars over the hemisphere. 














RADIATION PRESSURE IN GALACTIC NEBULAE 
JESSE L. GREENSTEIN 


ABSTRACT 
_ _ The existence of small solid particles with radii approaching the wave-length of light 
is indicated by the observed interstellar selective absorption. The recent investigations 
of Struve and his collaborators indicate that particles much larger than the wave-length 
of light dominate in bright galactic nebulae. The source of this difference is found in the 
enormous radiation pressure of the hot stars illuminating the nebulae. Figures 1-4 show 
the dependence on spectral type of the sizes of the particles repelled by a star. Further 
observations should yield information on the frequency function of particle sizes in 
space. 
I. DISTRIBUTION OF PARTICLE SIZES IN INTERSTELLAR SPACE 

The existence of selective absorption in the galactic plane indi- 
cates a considerable abundance of particles with radii of the order 
of, or less than, the wave-length of light. The fact that a low value 
is suggested by recent investigations” * for the ratio of the selective 
to the total absorption, and the spectrophotometric measurements 
of the distortion of the continuous radiation in early-type stars,* * 5 
prove that Rayleigh scattering cannot be an important process. 
Some observational evidence on the frequency of particles of various 
radii is available from the distribution of interstellar hyperbolic me- 
teors.° Over the observed range from 2.9 X 107" cm to 6.6 X 1073 
cm, Watson has shown that the frequency of meteoric particles in 
space increases uniformly with decreasing particle size.’ The total 
absorption produced by such particles may be considerable, but 
the reddening will be zero. Particles capable of producing selective 
absorption would be meteors fainter than the sixteenth magnitude. 


* Stebbins and Whitford, Ap. J., 84, 132, 1936. 

2 Greenstein, Harvard Ann., 105, 17, in press. 

3 Trumpler, Pub. A.S.P., 42, 267, 1930. 

4Struve, Keenan, and Hynek, Ap. J., 79, 1, 1934. 

5 J. Rudnick, ibid., 83, 394, 1936. 

6 Watson, Harvard Ann., 105, 32, in press. 

7 The frequency and the velocity distribution of hyperbolic meteors will undoubted- 
ly prove of great importance in the study of the interstellar medium. 


242 








































RADIATION PRESSURE IN GALACTIC NEBULAE 243 


The existence of selective absorption in the galaxy requires that the 
frequency function continue to increase at least as far as particles 
having radii of the order of the wave-length of light. The total 
density in space of these smaller particles may be less than that of 
the large particles, even though they dominate in the production of 
the total and selective absorption. The space density, p;, of par- 
ticles of radius a;, required to produce a given total absorption 
coefficient, is of the form: 


aj 


aie K(a;) ; 


K(a;) is an efficiency factor, nearly unity except for a; < (A/2m), 
when it is of the form (a/A).4 

The recent work of Struve and his collaborators,* ° Collins’ 
and other investigators on the colors of galactic nebulae with con- 
tinuous spectra has yielded some interesting information as to the 
distribution of particle sizes in the neighborhood of hot stars. The 
observed agreement of the colors of these nebulae with the colors of 
the illuminating stars shows that the dominant process is diffuse 
reflection, or nonselective scattering. Any considerable dependence 
of scattering on wave-length for any large fraction of the scattered 
light, such as would exist for a <X, is excluded by this observation. 
Struve and Story" have also shown that simple diffuse reflection 
from large particles (a > X) will satisfactorily explain the observed 
relation between the brightness of the central star and the angular 
diameter of the nebula, found by Hubble.” This remarkable devia- 
tion from the form of the frequency function suggested above for 
interstellar space, in luminous nebulae which are associated with 
apparently normal obscuring clouds, can be simply interpreted as 
a result of the enormous radiation pressure exerted on small par- 
ticles by hot stars of high luminosity. 


8 Struve, Elvey, and Keenan, Ap. J., 77, 274, 1933. 
9 Struve, Elvey, and Roach, ibid., 84, 219, 1936. 
10 Pub. Amer. Astr. Soc., 8, 240, 1936. 

1 Ap. J., 84, 203, 1936. 

12 Tbid., 56, 162 and 400, 1922. 








JESSE L. GREENSTEIN 


II. RADIATION PRESSURE 

The classical theory of the scattering of light by small, solid, 
spherical particles, developed by Mie," has been applied by Debye™ 
to the computation of the radiation pressure on such particles. 
Computations of the ratio of pressure to gravity, (P/G), by Gerasi- 
movic¢," as well as by Schoenberg and Jung,’ have shown that this 
ratio can reach very high values (~1o‘). The value of the pressure 
varies with the ratio of the size of the particle to the wave-length, 
and is given by 


wa?ky 27rd 
P= TF oa); a = 


Fy is the intensity; g(a) has been computed by Debye for several 
substances of differing optical properties, and depends on the index 
of refraction, v, and the absorption coefficient of the substance. 
Schoenberg and Jung have given a table of (P/G) for a perfect 
conductor (perfect reflector), vy = ©. Their table for radii less 
than 10~° cm gives too large values for the ratio (P/G) because of 
erroneous values they have adopted for g(a) for small values of a. 
Debye’s theory gives a simple expression for very small particles, 
which was used in the present computations: 


On the assumption that the star radiates as a black body, the 
value of (P/G) for a particle of radius a and density s is 


. 3Rz 7? 
= ‘ *K r j 
(6) re | E,¢(a)dr 


R,, and m,, are the radius and mass of the star; Gis the gravitational 

constant; c is the velocity of light; and £, is the black-body energy 

corresponding to the effective temperature, Te. The values of 

(P/G) have been computed by numerical integration in a wide 

range of particle sizes and stellar temperatures for the substances 
13 Ann. d. Phys., 25, 377, 1908. %Z.f. Ap., 4, 205, 1932. 


14 [bid., 30, 57, 1909. 16 A.N., 247, 413, 1932. 
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with the optical and physical properties given in Table 1. The 
errors arising from the lack of accurate values of g(a) and from the 
inaccuracies of the numerical integration are not larger than 10 per 
cent for most of the values of (P/G), except for v = 1.1. In any 











TABLE 1 
Refractive Inde: Density ; 
ee ee oi Possible Substance 
(v) (s) 

T / ; . 

(1—V —1) 5.0 Metals—iron, nickel, copper 
2 
oo4 0.5 Ice—Solid hydrogen, helium, or oxygen 
eee ree 2:0 Silica—stony meteors 
Infinity. 3.0 Perfect reflectors—‘‘white body” 











case, the results must be viewed as only approximate; the theory is 
developed for ideal smooth spherical particles with optical proper- 
ties independent of the wave-length. The deviation for very small, 
coarse aggregates such as meteoric dust may be serious. 








TABLE 2 
Effective | 
. or kr m, R, 
Spectral Type Temperature ps R 
(Degrees) © © 
je 25,000: 19.5: 7.8: 
Bi. =. 22,000: 13.0: 5.8: 
Base, 18,000: 8.6: 4.6: 
—_———e 15,500 6.2 3.9 
LS ae 13,000 4.7 4.4 
|. > ee etn 9400 2.6 2.4 
A5 MPL ee ICE 7700 1.9 1.8 
gs ; : 6000: 2.8: $8: 
gGs.. RTA 4500 06) 54.3: 
MONO xa we 3500 | 6.0: 63.9: 











In order to determine the ratio (P/G) for typical stars of the 
spectral sequence, it was necessary to use values for the effective 
temperatures, radii, and masses of the stars. The adopted values in 
Table 2 are taken from a recent compilation of these data by Pi- 
lowski."7 


17 Pilowski, A.N., 261, 18, 1936 
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The detailed values of (P/G) so obtained will be published and 
discussed elsewhere, together with a treatment of the general prob- 
lem of the stability of small particles in the galactic system. In the 
present discussion the results are given schematically for the various 
substances in Figures 1-4. The shaded areas in the figures indicate 
the range of particle radius, a, within which the ratio of pressure to 
gravity, (P/G), is greater than unity, as a function of spectral type. 
Particles for which this ratio is greater than unity cannot move in 
closed orbits around the central star, and will therefore be infre- 
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quent, if not entirely absent, at distances that are large compared 
with the usual sizes of luminous nebulae. 

The ratio (P/G) depends on the total luminosity of the star and 
also on the wave-length of the energy maximum. The decreasing 
luminosity along the main sequence is the obvious origin of the fea- 
tures of that part of the diagrams. The late-type giants, though of 
high luminosity, have energy maxima strongly displaced toward 
the longer wave-lengths. The value of g(a) for a given particle size 
is therefore smaller for the red stars, reducing the pressure sharply 
for small particles. For metallic particles, the ratio (P/G) is inde- 
pendent of the radius for all particles of radius less than about 
10-> cm. 











RADIATION PRESSURE IN GALACTIC NEBULAE 


III. EFFECTS OF RADIATION PRESSURE IN GALACTIC NEBULAE 

From these considerations it is obvious that the frequency func- 
tion of particle sizes in space will be seriously affected near the very 
stars which most commonly excite nebulae. Particles whose radii 
nearly correspond to the wave-length of light are most subject to this 
effect and are absent from the nebulae. An examination of the figures 
indicates two general effects: (1) The total space density of inter- 
stellar matter will be decreased in the neighborhood of hot stars to a 
degree which depends on the spectral type. (2) Those particles 
capable of producing selective absorption (log a ~ —4.6) are 
excluded from the neighborhood of stars of types Bo to As, and gF5 
to gMo."* Very small particles (log a < —6.0) may be present 
near stars of types Ao and later, and of types gF5 to gMo. 

1. Struve and Story’® have shown that the total obscuration by 
dark material near emission-line nebulae is in fact less than that 
near reflection nebulae, indicating that the space density is lower 
near the hotter stars, in accordance with the theoretical predictions. 
The density 01 the scattering material present near emission nebulae 
must be less than that near reflection nebulae in order to explain 
the fact that the surface brightnesses of emission-line nebulae were 
no greater than those of reflection nebulae, in spite of the large 
amount of additional radiation in the strong nebular emission lines. 
This point was first noted by Struve. 

2. The nonappearance of selective scattering in the nebulae so 
far observed may be due to the complete absence of small particles 
from the dark clouds with which the bright nebulae are associated. 
This is rendered rather improbable by the observed reddening of 
objects viewed through dark nebulae. The theory developed above 
will account simply for the distortion of the frequency function that 
is indicated by the agreement of the color of the nebulae with that 
of the illuminating star. A further prediction subject to observa- 
tional test is that nebulae associated with stars as late as Ao, or 
with late-type giants, might show selective scattering and be 

18 Gerasimovit, Ap. J., 78, 298, 1933, has made a similar suggestion on the basis of 


less complete calculations. 


t9 Loc. cit. 
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bluer than the illuminating star. The observational results are not 
definite as yet; the nebula near Antares? (cMo) seems to be nearly 
of the same color as Antares, while the nebula near y Cygni? (F8p) 
is bluer than the star.”? Such small particles are very inefficient 
scattering agents per unit mass, since the efficiency factor K(a;) is 
of the form (a/X)*. Therefore, only if the space density is very high 
will the nebula show strong selective scattering. 

The importance of observations of the colors of nebulae near 
late-type stars cannot be overemphasized. Since space reddening 
may seriously affect some of the measured colors of nebulae, it is 
advisable to determine the color indices of the illuminating stars 
in the same manner as that of the nebulae. The frequency function 
for small particles can be investigated in this manner far beyond the 
range of sizes observable as meteors. Such data will prove of great 
importance in the discussion of the problems connected with the 
galactic absorption coefficient and the physical nature of the inter- 
stellar medium. 


HARVARD COLLEGE OBSERVATORY 
February 12, 1937 


20 A private communication from Struve and Elvey indicates that they suspect the 
existence of Ha emission in the nebula near y Cygni. Therefore the relative blueness 
of the nebula cannot be simply interpreted since other emission phenomena may exist. 





































NOTES 


LINE STRENGTHS IN NEON I 
ABSTRACT 

The line strengths for the transition array 2p54p— 2p54s of neon I have been calcu- 
lated in intermediate coupling, including the electrostatic and spin-orbit interactions. 

The line strengths for the transition array 2p53p— 2p*3s of neon I 
have been calculated by G. H. Shortley' in intermediate coupling, 
using the electrostatic and spin-orbit interaction terms and evaluat- 
ing the radial integral parameters from the observed energy values. 
Similarly, the results of calculating the line strengths for the neon 
transition array 2p°4p— 2p*4s are given here. 

Neon 2p‘4s fits Houston’s formulas well, so that the error in the 
line strengths introduced by the parameters from this configuration 
is probably less than 1 per cent. The radial integral parameters for 
neon 2p'4p have been calculated by C. L. Bartberger? in obtaining 
the energy-levels. Shortley found that a small change in the values 
of these parameters made a large change in the line strengths. How- 
ever, as the calculated values of the energy-levels for 2p°4p agree 
better with the observed values than those for 2p*3p, the radial 
integral parameters for 2p*4p should give more accurate values of 
the line strengths. The greatest deviation of a calculated level from 
an observed one was 70 cm~ for 2p*3p, 30 cm™ for 2p°4p; and the 
sum of the squares of all the deviations was 13,900 for 2p°3p, 2900 
for 2p54p, where 2p53p has a spread of 4713 cm between its lowest 
and highest energy-levels and 2p%4p a spread of 1768 cm™. It is 
convenient to estimate the error in calculating energy-levels by the 
ratio of the sum of the squares of the deviations to the spread of the 
configuration, which gives 294 per cent for 2p°3p and 164 per cent for 
2p*4p. 

The values of the line strengths S(p, s) are given in Table 1. 
These are not the experimental intensities’ of the lines but are re- 

t Phys. Rev., 47, 295, 1935. 2 [bid., 48, 682, 1935. 

3 E. U. Condon and G. H. Shortley, The Theory of Atomic Spectra (Cambridge Uni- 
versity Press, 1935), § 7, pp. 97 and 240. 
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lated to the spontaneous transition probabilities A(p, s) of Einstein 
by the formula S(p, s) = 3hc3(2J, + 1) A(p, s)/647‘v3. vy is the 
frequency of the line, and the values of J, are given as subscripts on 
the right of the Russell-Saunders notation for 2p54p. The Russell- 
Saunders level which occurs in largest proportion in each of the levels 
is used to denote that level. In addition, the original notation of 
Paschen and G6tze is given for reference. 


TABLE 1 


LINE STRENGTHS IN THE TRANSITION ARRAY 2p54p— 2p54s, NEON I 
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C. W. Urrorp 
ALLEGHENY COLLEGE 
February 1937 


NOTE ON PHOTOELECTRIC PHOTOMETRY 

In a recent paper’ I published observed magnitude differences be- 
tween Vega and 8 Pegasi as a function of the mean wave-length. 
This mean wave-length corresponded to that at the center of the 
slits and not to the effective wave-length. Since the slits admitted 
480 A at a time, the difference between the mean and the effective 
wave-length is not always negligible. 

Smooth curves were drawn through the observed points in Figure 
1.2 The abscissa of this graph is the mean wave-length, and the ordi- 
nate is the electrometer deflection or the effective energy. By a 
simple graphical process the effective wave-length was found for each 


1 Ap. J., 84, 369, 1936. 2 Tbid. 
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star for each mean wave-length. Since the shape of each curve is 
well defined, it is possible to correct the energy observed at the effec- 
tive wave-length to the value it would have at the mean wave- 
length. This process was carried through for both stars. 

It was stated in the former paper* that these two radically differ- 
ent stars were compared primarily to test the efficacy of the appa- 
ratus in the determination of the effective wave-lengths when two 
color filters were used with the apparatus. Although these effective 
wave-lengths are not substantially changed by the corrections men- 
tioned above, it seems worth while to give the corrected magnitude 
differences between these two stars. 

In the following table the mean magnitude differences given in 
Table I¢ have been corrected in such a way that the wave-lengths 
are now effective wave-lengths. More than half of the corrections 
were 0.02 mag. or less, and none was more than o.10 mag. 








| EFFECTIVE WAVE-LENGTH (yu) 
©.454 0.478 0.501 0.525 0.549 0.573 0.597 0.621 
6B Peg-a Lyr +363 | 3Moo 2Mgo 2Mc4 230 2Mog 1Mo2 169 
CS) eee ....{ 0.645] 0.669 | 0.693 | 0.717 | 0.741 | 0.765 | 0.789 0.813 
AM. cecal ExQ@@h Bo2i 1.07 °.79 0.47 0.31 0.18 |+0.06 
(u) eee 0.838] 0.862 | 0.886 | 0.910 | 0.934 | 0.958 | 0.983 I .007 
I 2 eee ee —0.06 | 0.15 0.17 0.30 0.23 0.36 0.44 |—0.52 





























As a result of these corrections the observed difference in the color 
temperatures of these two stars is slightly larger. It still appears 
that if 8 Pegasi has a black-body temperature of 3000° K, the cor- 
responding temperature of Vega is between 10,000° and 13,000", de- 
pending on the spectral regions compared. If 6 Pegasi is hotter than 
3000°, Vega is hotter than the foregoing figures indicate. This last 
statement was erroneously given in the original paper. 

Joun S. HALL 


SPROUL OBSERVATORY 
SWARTHMORE COLLEGE 
SWARTHMORE, PENNSYLVANIA 


3 bid. 4 Tbid. 
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OBSERVATIONS OF THE SPECTRUM OF 
NOVA HERCULIS 1934 

At the request of the editors of the Astrophysical Journal the 
directors of a number of observatories kindly communicated de- 
tails of the spectrograms of Nova Herculis 1934 secured at their 
observatories between December 13, 1934, and April, 1935. On the 
formation of a subcommission on novae in the Commission on 
Stellar Spectra of the International Astronomical Union these de- 
tails were passed on to the subcommission. They have been recently 
completed by returns from some outstanding observatories. Alto- 
gether 34 observatories obtained well over 3000 spectrograms be- 
tween December 13, 1934, and April 20, 1935; and on every day in 
that interval at least one spectrogram was secured at some observa- 
tory. Most observatories are working up their own material, but 
in neverly every case the directors have expressed their willingness 
to send copies of plates to those investigators who are anxious to 
clear up difficulties arising from gaps in their own series. The com- 
plete list of observations is too long to publish, and a summarized 
statement may be more misleading than helpful. But the data are 
in my possession; and if any investigator wishes to know where he 
can find material to fill in a gap in his series of plates, I will, on 
receipt of a letter stating the dates and the region of the spectrum 
desired, send all relevant information in my possession. 

F. J. M. STRATTON 
GONVILLE AND CAIUS COLLEGE 
CAMBRIDGE, ENGLAND 


COLORS OF NEBULAE NEAR y CYGNI 

In a recent paper’ the color index of the diffuse nebula near 
BD+39°4206 was found to be —o.2 mag. If this nebula is produced 
by nonselective scattering of the light of y Cygni, the color index 
should be roughly +o0.5 mag. Since Henyey? found no measurable 
polarization in the light of this nebula, Rayleigh scattering by very 
small particles cannot account for its apparent blueness. In our de- 
terminations of the color index we used Eastman emulsion IG with 

tAp. J., 84., 226, 1936. 

2 Ibid., p. 617. 
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a yellow filter (Wratten No. 12, “minus blue’). This filter cuts off 
quite sharply at A 5000, while emulsion IG cuts off at approximately 
5800. This spectral region does not include any of the hydrogen 
lines, and the nebular line 5007 lies in the region where the filter 
begins to absorb appreciably. An emission nebula would appear to 
be quite faint if photographed with this combination. For the violet 
region we used Eastman Process or lantern-slide plates, which cover 
the region from \ 3000 to about A 4700. All hydrogen lines from the 
Balmer limit to Hy and many strong nebular lines are therefore in- 
cluded. In order to test whether the nebulae near y Cygni are of the 
emission type we have taken several photographs on emulsion IC 
with a red filter (Wratten No. 23A). This combination includes the 
line Ha. The three nebulae near BD+41°3687, BD+39°4206, and 
BD+ 39 4186 are quite strong in red light. The spectral type of 
BD+41°3687 is Ko, but since no conspicuously red stars are near 
the other two nebulae, we tentatively infer that all three nebulae 
are of the emission type. If, as is probable, all three nebulae derive 
their light from the F8p star y Cygni, we are confronted with the 
interesting problem of accounting for the emission. Unfortunately, 
these nebulae are so faint that we have not been able to obtain 
their spectra with the objective prism. 
OTTO STRUVE 
C. T. ELVEY 
March 1937 


THE MAGNITUDE OF THE COMPANION 
TO PROCYON 

The companion to Procyon resembles in many ways the com- 
panion of Sirius: (1) the history of the discovery is similar; (2) the 
companions are very faint, compared with their primaries (we find 
that Am = to mag. in both cases); (3) both companions are highly 
underluminous, as compared with ordinary stars of the same mass. 
This last feature indicates that both stars are white dwarfs, which 
is verified for the Sirius companion by its spectrum. 

From the experience with gratings in measuring Am in close 
double stars' the writer was convinced that the best method of 
t Pub. A.S.P., 46, 99, 1934; 47, 24 and 166, 1935 
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measuring the brightness of the Procyon companion would be a de- 
termination of Am with the aid of a double grating.’ 

At first a grating was made with a free opening of about 34 inches 
square, which gave first-order images 5.9 mag. fainter than the 
central image. When placed on the 40-inch telescope on a good night, 
this grating appeared to produce (1, 
1) images too faint to be seen (Am 
B- @* = 2X 5.9 = 11.8 mag.), whereas the 
companion was visible. Accordingly a 








new grating was made, giving Am = 2 X 


Xi A 5.1 = 10.2 mag. The wires are 0.1150 

e © @ inches in diameter and are spaced by 
Saal - ae 1.300 inches (center to center), giving 
. Am = 10.186 mag. Twenty-five wires 
= = are placed in each co-ordinate, giving 








a square opening of 26 X 1.3 = 33.8 
inches. The visual appearance of the 
patterns without and with grating is 
shown in Figure 1. Comparisons were made, with and without a 
yellow filter, between the images a and B, giving 





FIG. 1 


B=a+to™r +02, 
Am = 10.3 + 0.2 (m.e.). 


The mean error is an estimate; it is considered very improbable that 
Am lies outside the limits 10.0 mag. and 10.5 mag. 

The apparent magnitude of Procyon B is therefore 10.8, and the 
absolute magnitude is 13.1, not nearly as faint as was often sup- 
posed. 

The measurement of Procyon completes the writer’s Am deter- 
minations of the (roughly five hundred) visual binaries brighter 
than 6.50 mag., closer than 5"0, north of —2o0° declination. The 
extension of the distance limit to 30” is being undertaken. 

G. P. KUIPER 
YERKES OBSERVATORY 
March 29, 1937 

2It turned out that Am was considerably smaller than expected, and there would 

have been no objection to using a single grating; no very thin wires would have been 


needed. 
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A NEW BINARY OF LARGE PARALLAX 


In the systematic micrometer work on the stars of large parallax, 
BD —4°3665 was found to have a distant companion which shares 
its proper motion, 0766 in 261°. This is shown by the following 


measures: 
i) a? ee 107°3 15°34 m(B)=14.3 36, 540X seeing 3+ 
SONG. 68. . 108.3 15.23 15.3 36, 360 2+* 
1937.27...... 107.0 15.35 14.8 40, 350 3 
£046). GE. <0 0: 107.4 35.32 14.8 3n 
* Half-weight, faint. 
The position for 1900 is a= 14514™4,6= — 4°41’. The magnitude of 


A is 7.51 on the international photovisual scale, as determined from 
two sets of measures made at Leiden. The Mount Wilson spectrum 
is dK1r, and the spectroscopically determined absolute magnitude 
is 5.5. The spectroscopic parallax is therefore o%o40. 

The trigonometric parallax is 07053 + .006 according to Schles- 
inger. Adopting this value, the absolute visual magnitude of the 
companion becomes 13.4; the spectroscopic parallax would make it 
12.8 mag. The companion is, therefore, one of the faintest stars 


known. 
G. P. KUIPER 
YERKES OBSERVATORY 
April 14, 1937 


NOTE ON INTERSTELLAR SCATTERING 

In a recent paper by Elvey and Roach’ on the light from the 
night sky it was shown that the intensity of the light scattered by 
interstellar matter is very nearly the same as the intensity of the 
light coming to the earth from all the stars. In analyzing their 
data these writers made use of an approximate formula, derived by 
Struve,? giving the relationship between starlight and scattered 
light. It is the purpose of this note to derive an expression which 


is more rigorous. 


1 Ap. J., 85, 213, 1937. 2 Tbid., p. 194 
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The equation of transfer governing the problem is 


dI 


dr 


+ /] = [Ibde + { Tobde : 


where J and J, are the intensities of the scattered light and the 
starlight, respectively. Now if I, is taken to be everywhere equal 
to L, where L is the average intensity of starlight as seen from the 
earth, and if J is assumed isotropic, we have 


dl 


F a ee, 


where ¥ is written for the albedo {@dw. The solution of this differ- 


ential equation is 


When + => © and when y is small, this reduces, to within terms in 
y’, to the approximate expression referred to. 
Since J and L are the same, the condition imposed on y and 7 


by the observations is 


Table 1 gives the values of 7 for the values of y compatible with 


the observations. 


TABLE 1 
Y T ¥° T 
OU. Sos co mag QO. aes TE mae 
OO es ei 3.0 O-0:. E32 
OF icles 2.0 FRO so es Pal 


It will be seen that the albedo must, within the errors of observa- 
tion, not exceed 0.5. Since it is not sensitive to the optical depth 
when the latter is large, it is likely that y is about o.5 or 0.6. 

L. G. HENYEY 
YERKES OBSERVATORY 
WILLIAMS Bay, WIs. 
April 1937 


